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ABSTRACT 



The scope of this thesis is to evaluate the use of the MINI RANGER, 
LORAN (TD) and LORAN (DISPLAY) navigation systems in order to support 
the collection of current profiles by an Acoustic Doppler Current Profiler Re- 
corder (ADCP). 

A theoretical error analysis of these systems is undertaken in order to estab- 
lish the minimum error limits which might be expected when averaging current 
profiles over time frames of up to 30 minutes. h , 

Experimental data was collected with all of these systems in the Monterey 
Bay and was analysed, the results being presented in this thesis. In addition GPS 
data was also collected but time has prohibited its analysis and subsequent 
inclusion. 

The results show that because of ship fluctuations in course and speed there 
is no statistically significant difference between the navigation systems when we 
average the data o\'er times of about 25 - 30 minutes. However, they also show 
that both the MINI RANGER and Bottom Tracking with the ADCP can 
produce reasonable results in as little as three minutes, although the ADCP re- 
sults are clearly biased. 



TABLE OF CONTENTS 



I. INTRODUCTION 1 

II. DATA COLLECTION 2 

A. CRUISE 2 

B. TIME SYNCHRONIZATION 2 

III. ACOUSTIC DOPPLER CURRENT PROFILERS 5 

A. INTRODUCTION 5 

B. PRINCIPLES OF OPERATION OF THE ACOUSTIC DOPPLER 

CURRENT PROFILER 6 

C. THREE DIMENSIONAL CURRENT VELOCITY VECTORS ...8 

D. DEPLOYMENT AREA 9 

E. BOTTOM TRACKING OPERATION 10 

F. VELOCITY PROFILES 11 

G. ACOUSTIC DOPPLER CURRENT PROFILER DATA PROC- 
ESSING 11 

IV. THEORY OF LORAN C 13 

A. INTRODUCTION 13 

B. SYSTEM DESCRIPTION 13 

C. PULSES-PHASE AND CODES-CYCLE SELECTION 13 

D. TIME DIFFERENCE EQUATIONS 16 

E. HYPERBOLIC GRADIENT, CROSSING ANGLES 18 

F. ERROR ANALYSIS FOR LORAN C 20 

V. THEORY OF MINI RANGER 30 

A. INTRODUCTION 30 

B. REFERENCE STATIONS 30 



C. RANGE POSITIONING 30 

D. ERROR ANALYSIS FOR MINI RANGER 31 

VI. RESULTS 36 

VII. CONCLUSIONS AND RECOMMENDATIONS 46 

APPENDIX A. STATISTICAL ANALYSIS FOR THE LC 408 DATA 47 

APPENDIX B. TABLES FOR STANDARD DEVIATIONS AND 
COVARIANCES OF VELOCITY COMPONENTS FOR LORAN C (LC 
408) 48 

APPENDIX C. PROGRAM DRIVER FORTRAN 51 

A. GENERAL REMARKS 51 

B. DESCRIPTION OF PARAMETERS; 51 

APPENDIX D. TABLE FOR STANDARD DEVIATIONS AND 
COVARIANCES OF VELOCITY COMPONENTS FOR MINI RANGER 66 

APPENDIX E. AVERAGE CURRENT PROFILES 67 

APPENDIX F. CURRENT PROFILES FOR SPECIFIC TIMES 79 

LIST OF REFERENCES 85 



INITIAL DISTRIBUTION LIST 



86 



LIST OF TABLES 



Table 1. COMPARISON OF CLOCKS WITH THE MINI RANGER 

(MR) 2 

Table 2. DIFFERENT TIMES BETWEEN THE CLOCKS WITH COM- 
MON TIME BASE (ADCP) 2 

Table 3. LORAN C PHASE CODES 13 

Table 4. REFERENCE STATION SITES 30 

Table 5. POSITION ERROR FOR DIFFERENT CROSS ANGLES . . 32 

Table 6. MAX-MIN VALUES OF STANDARD DEVIATIONS AND 
COVARIANCES FOR VELOCITY COMPONENTS FOR 3 
MINUTE TIME INTERVALS BETWEEN POSITIONS (MINI 

RANGER) 34 

Table 7. SHIP VELOCITY COMPONENTS FOR MINI RANGER, 

LORAN (TD), LORAN (DISPLAY) FROM POINT 1 TO POINT 

2 AT 3 MINUTE TIME INTERVAL 36 

Table 8. SHIP VELOCITY COMPONENTS FOR MINI RANGER, 

LORAN (TD), LORAN (DISPLAY) FROM POINT 2 TO POINT 

1 AT 3 MINUTE TIME INTERVAL 36 

Table 9. COMPARISON BETWEEN MINI RANGER AND LORAN 
(TD) AND THE MINI RANGER AND LORAN (DISPLAY) 

FOR LEGS 1 - 2 AND 2-1 37 

Table 10. SHIP VELOCITY COMPONENTS FOR MINI RANGER, 

LORAN (TD), LORAN (DISPLAY) FROM POINT 1 TO POINT 

3 AT 3 MINUTE TIME INTERVAL 38 

Table 11. COMPARISON BETWEEN MINI RANGER AND LORAN 

(TD) AND THE MINI RANGER AND LORAN (DISPLAY) 

FOR LEG 1-3 39 

Table 12. SHIP VELOCITY COMPONENTS FOR MINI RANGER, 
LORAN (TD). LORAN (DISPLAY), BOTTOM TRACKING 



VI 



FROM POINT 3 TO POINT 2 AT 3 MINUTE TIME INTER- 
VAL 

Table 13. SHIP VELOCITY COMPONENTS FOR MINI RANGER, 
LORAN (TD), LORAN (DISPLAY), BOTTOM TRACKING 
FROM POINT 2 TO POINT 3 AT 3 MINUTE TIME INTER- 
VAL 

Table 14. COMPARISON BETWEEN MINI RANGER AND LORAN 
(TD) AND THE MINI RANGER AND LORAN (DISPLAY) 

FOR LEGS 3 - 2 AND 2-3 

Table 15. SHIP VELOCITY COMPONENTS FOR MINI RANGER, 
LORAN (TD), LORAN (DISPLAY), BOTTOM TRACKING 
FROM POINT 3 TO POINT 4 AT 3 MINUTE TIME INTER- 
VAL 

Table 16. COMPARISON BETWEEN MINI RANGER AND LORAN 
(TD) AND THE MINI RANGER AND LORAN (DISPLAY) 

FOR LEG 3-4 

Table 17. AVERAGE VELOCITY COMPONENTS FOR MINI 

RANGER, LORAN (TD). LORAN (DISPLAY), BOTTOM 

TRACKING 

Table 18. TIME INTERVAL BETWEEN POSITIONS 3 MINUTES . 
Table 19. TIME INTERVAL BETWEEN POSITIONS 6 MINUTES . 
Table 20. TIME INTERVAL BETWEEN POSITIONS 12 MINUTES 
Table 21. TIME INTERVAL BETWEEN POSITIONS 20 MINUTES 
Table 22. TIME INTERVAL BETWEEN POSITIONS 40 MINUTES 



39 



40 



40 



41 



43 



44 



48 

48 

49 

49 

50 



vii 



LIST OF FIGURES 



Figure 1. CHART OF THE WORK AREA 

Figure 2. VOLUME SCATTERING OF SOUND 

Figure 3. FIRST-SECOND SOUND DOPPLER SHIFT 

Figure 4. WATER VELOCITY VECTOR 

Figure 5. VELOCITY RESOLUTION WITH FOUR ADCP BEAMS . . . 

Figure 6. SMALL AND LARGE ERROR VELOCITY 

Figure 7. ADCP AND CURRENT METERS 

Figure 8. LOCATION OF WEST COAST LORAN C STATIONS 

Figure 9. LORAN C GROUP REPETITION INTERVAL (GRI) AND 

TIME DIFFERENCES (TD) 

Figure 10. THIRD CYCLE TRACKING POINT 

Figure II. TYPICAL HYPERBOLIC FIX GEOMETRY 

Figure 12. HYPERBOLIC GRADIENTS 

Figure 13. HYPERBOLIC GRADIENTS AND CROSSING ANGLES . 
Figure 14. SPHERICAL TRIANGLE SPANAGEL-M-NORTH POLE . 

Figure 15. PLANNING THE MOST DESIRABLE SITE 

Figure 16. MAXIMUM POSITION ERROR 

Figure 17. AVERAGE CURRENT PROFILE FROM POINT 1 TO 
POINT 3 USING NAVIGATION DATA FROM MINI 

RANGER 

Figure 18. AVERAGE CURRENT PROFILE FROM POINT 3 TO 
POINT 2 USING NAVIGATION DATA FROM MINI 

RANGER 

Figure 19. AVERAGE CURRENT PROFILE FROM POINT 2 TO 
POINT 3 USING NAVIGATION DATA FROM MINI 

RANGER 

Figure 20. AVERAGE CURRENT PROFILE FROM POINT 3 TO 
POINT 4 USING NAVIGATION DATA FROM MINI 



3 

6 

7 

8 

8 

9 

11 

14 

15 

16 

17 

19 

20 

21 

31 

33 

67 

68 

69 



RANGER 70 

Figure 21. AVERAGE CURRENT PROFILE FROM POINT 1 TO 

POINT 3 USING NAVIGATION DATA FROM LORAN (TD) 71 
Figure 22. AVERAGE CURRENT PROFILE FROM POINT 3 TO 

POINT 2 USING NAVIGATION DATA FROM LORAN (TD) 72 
Figure 23. AVERAGE CURRENT PROFILE FROM POINT 2 TO 

POINT 3 USING NAVIGATION DATA FROM LORAN (TD) 73 
Figure 24. AVERAGE CURRENT PROFILE FROM POINT 3 TO 

POINT 4 USING NAVIGATION DATA FROM LORAN (TD) 74 



Figure 25. AVERAGE CURRENT PROFILE FROM POINT 1 TO 
POINT 3 USING NAVIGATION DATA FROM LORAN 

(DISPLAY) 75 

Figure 26. AVERAGE CURRENT PROFILE FROM POINT 3 TO 
POINT 2 USING NAVIGATION DATA FROM LORAN 

(DISPLAY) 76 

Figure 27. AVERAGE CURRENT PROFILE FROM POINT 2 TO 
POINT 3 USING NAVIGATION DATA FROM LORAN 

(DISPLAY) 77 

Figure 28. AVERAGE CURRENT PROFILE FROM POINT 3 TO 
POINT 4 USING NAVIGATION DATA FROM LORAN 

(DISPLAY) 78 

Figure 29. CURRENT PROFILE FROM POINT 1 TO POINT 2 USING 
NAVIGATION DATA FROM MINI RANGER FOR A SPE- 
CIFIC TIME (14 33 II) 79 

Figure 30. CURRENT PROFILE FROM POINT 2 TO POINT 1 USING 
NAVIGATION DATA FROM MINI RANGER FOR A SPE- 
CIFIC TIME (15 06 10) 80 

Figure 31. CURRENT PROFILE FROM POINT 1 TO POINT 2 USING 
NAVIGATION DATA FROM LORAN (TD) FOR A SPE- 
CIFIC TIME (14 33 1 1) 81 

Figure 32. CURRENT PROFILE FROM POINT 2 TO POINT 1 USING 



IX 



NAVIGATION DATA FROM LORAN (TD) FOR A SPE- 
CIFIC TIME (15 06 10) 

Figure 33. CURRENT PROFILE FROM POINT 1 TO POINT 2 USING 
NAVIGATION DATA FROM LORAN (DISPLAY) FOR A 

SPECIFIC TIME (14 33 11) 

Figure 34. CURRENT PROFILE FROM POINT 2 TO POINT 1 USING 
NAVIGATION DATA FROM LORAN (DISPLAY) FOR A 
SPECIFIC TIME (15 06 10) 



I. INTRODUCTION 



The relative accuracy of navigation data is very important in the processing 
of Acoustic Doppler Current Profiler (ADCP) data. 

The shipboard ADCP is an instrument which measures vertical profiles of 

r 

currents. The raw ADCP data is collected in 128 four meter bins, and is used in 
conjuction with a navigation system from which the position of the ship is re- 
corded. 

In this experiment we will try to evaluate which of the navigation systems 
a\ailable to us is better for use with the ADCP in order to determine the absolute 
current velocity. We collected data with MINI RANGER, LORAN (TD), 
LORAN (DISPLAY), and GPS receivers. Due to time constraints, however, the 
GPS data has not been included in the analysis. 

Chapter 2 describes the cruise from which the data analysed in this thesis is 
taken. It also describes the time synchonization problems which existed between 
the various types of na\ igational equipment used. 

Chapter 3 describes the ADCP and its principles of operation. 

In chapters 4 and 5 we try to define the theoretical random error for the 
LORAN (LC 408) and the MINI RANGER respectively. 

Chapter 6 gi\es the results from the data collected in a cruise of Monterey 
Bay while Chapter 7 gives a summary of the conclusions from this study together 
with recommendations for future work. 
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II. DATA COLLECTION 



A. CRUISE 

The day chosen for the cruise was 22 September 1989 with the ship R/V Point 
Sur. During the cruise the following equipment belonging to the NPGS and R/V 
Point Sur was used; 

1. Motorola Mini Ranger, 

2. LORAN C receiver (Internav LC 408), 

3. Trimble GPS receiver, 

4. HP Computer logging three data sources simultaneously, and 

5. RDl ADCP with the Data Acquisition System (DAS) and IBM XT com- 
patible computer. 

The results reported in later sections will concern only the MINI RANGER 
and LORAN C navigational systems. 

The ship navigated the runs from point 1 to point 2, 2 - 1, 1 - 3, 3 - 2, 2 - 3, 
and 3 - 4 (see Figure 1). The speed during the measurements was 6 knots which, 
from pre\ious experience, has proven effective for ADCP data collection. In 
running from 1 - 2 for example, a turn was made at point 2 and then the ship 
steamed back to point 1 without attempting to follow the identical outward track. 

The a\ erage depth in the area of the triangle 1 2 3 is 30 fathoms, and for the 
run 3 - 4 it varies between 20 to 300 fathoms. 

B. TIME SYNCHRONIZATION 

During the cruise the four pieces of equipment used each had clocks which 
operated independently of each other. In order to sjmchronize all these clocks, at 
five times during the cruise time comparisons were made between them by using 
the Mini Ranger as the common time base (see for example Table 1). 

These times were then converted from a MR time base to an ADCP time 
base (see Table 2). 

The standard error in the time comparisons was at least 1 second. Once a 
dedicated multi-channel data logger is installed which allows all equipment out- 




Figure 1. CHART OF THE WORK AREA 

puts to be recorded on a common time base, it is anticipated that this synchroni- 
zation error will be reduced to negligible levels. 
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Table 1. COMPARISON OF CLOCKS WITH THE MINI RANGER (MR) 



HP 


15:16:28 


MR 


15:16:28.6 


LC 408 


15:14:00 


MR 


15:14:00.6 


ADCP 


15:15:10 


MR 


15:15:06.6 



Table 2. DIFFERENT TIMES BETT^ EEN THE CLOCKS WITH COMMON 



TIME BASE (ADCP) 



HP 


15:16:28 


ADCP 


15:16:32.0 


LC 408 


15:14:00 


ADCP 


15:14:04.0 


MR 


15:15:06.6 


ADCP 


15:15:10 



With five tables on a common ADCP time base, a linear regression between 
the different times was performed for each equipment type in order to calculate 
the coefficients A and B of the equation: 



Y=AX+B (1) 

The linear clock drifts were thus calculated and used to assist in interpolating 
the data into the ADCP time base. 



4 



III. ACOUSTIC DOPPLER CURRENT PROFILERS 



A. INTRODUCTION 

The shipboard Acoustic Doppler Current Profiler (ADCP) from RD Instru- 
ments is a state of the art instrument which measures vertical profiles of currents 
of the ocean at points along the ship's path [RD Instruments, 1989]. It allows the 
ocean to be sampled in a way which is fundamentally different from moored 
currents or drifters (instruments most commonly used for current measurement). 

Consider the general case of measuring the current at the water parcel whose 
coordinates are using an instrument whose location is A'o(0 



= ^^o(0 + ^ 



Where r is a vector between the two locations and 



-■ = a ;,(0 - A ^(/) 



From (2) and (3) 



a;,(/) = A’o(/) + [A'JO - Ao(0] 



The current U at the measurement point Ag -f- r will be: 



L’(A’o 4- /•) = 




<^-^0 ^(An- — Aq) 

dt ^ dt 



( 2 ) 

(3) 

(-0 

( 5 ) 



Equation (5) can be rewritten as 

U(X„ + r) = ^+ V(r) (6) 

Where V(r) is the velocity of the water parcel relatit e to the instrument posi- 
tion and dXo dt is the \ clocity of the instrument itself with respect to the Earth. 

Since the currents are a difference between two directly measured quantities 
dXo dt and V. the measurement of currents to an accuracy of 1 cm sec from a 



ship travelling at 10 knots (approximately 5.2 m/'see or 520 eni/see) requires that 
both dXo/dt and V must be measured to an aeeuraey of 0.2%. 

B. PRINCIPLES OF OPERATION OF THE ACOUSTIC DOPPLER 
CURRENT PROFILER 

ADCP's use the Doppler effect by transmitting a succession of aeoustie pulses 
at a fixed frequeney and listening to the resulting baekseattered water mass ech- 
oes in as many as 128 depth cells (bins) over a depth range of 30 to 700 meters. 
When seatterers move toward the ADCP, the echoes heard by the seatterers is 
Doppler shifted to a higher frequency. The amount of this shift is proportional 
to the relative veloeity between the ADCP (ship) and seatterer Figure 2. 




Part of these Doppler-shifted echoes refleet backwards or are baekseattered 
to the ADCP. The baekseattered echoes appear to the ADCP as if the seatterers 
were the souree of the eehoes (see Figure 3). 

The ADCP hears the baekseattered echoes, Doppler shifted a second time. 
Since the ADCP both transmits and receives, the Doppler shift is doubled, and 
the equation for the Doppler shift will be 

Fo = 2Fs-^cosA (7) 



where 

• Fj) is the Ddppler shift. 
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• Fs is the frequency of the sound when everything is still, 

• A is the angle between the acoustic beam and the water velocity, 

• V is the relative velocity between the sound source and the sound receitci, 
and 

• C is the speed of sound (m/sec) in the water at the face of the transducer, 
and can be calculated from the expression 

C = 1449.2 + 4.67- 5.5 lO'^r^ + 2.910-"7^ 

+ (1.34 - 10“^70(S - 35) + 1.610"^I> 



where 

• D = Depth, in meters, 

• S = Salinity, in practical salinity units, and 

• T = Temperature, in degrees Celsius. 

The water temperature at the face of the transducer is measured by the 
ADCP. The speed of sound (C) can be calculated from this measured water 
temperature using nominal values for depth (D) and salinity (S). 
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Computer analysis of the Doppler frequency shift of backscattered echoes 
from each bin is used to generate a precise depth segmented picture (or profile) 
of water currents throughout the water column bounded by the path of the 
acoustic beams. 

C. THREE DIMENSIONAL CURRENT VELOCITY VECTORS 

The ADCP beams each measure a single velocity component, i.e. the compo- 
nent of velocity toward or away from the transduser. When using multiple 
beams, one must make an assumption that currents are the same (homogeneous) 
over layers of constant depth. When the ADCP uses multiple beams pointed in 
different directions, it senses different velocity components. 




8 



Figure 5 shows how the ADCP, using four acoustic beams, computes three 
velocity components. The first pair of beams produces one horizontal and one 
vertical component, the second pair produces one horizontal (perpendicular to the 
first horizontal) and again one vertical component. 

The product SCALE y. {V^ + — Fj — F 4 ) is defined as the error velocity 

component, where F, to V^ are the velocity components along the beam directions 
Figure 6 . This sum should be close to zero. The error velocity component allows 
us to evaluate if the assumption for the horizontal homogenuity is reasonable. 
Also this is another way to estimate data quality. 




In Figure 6 there are two different situations, in the first the velocity is the 
same in all four beams at a constant level, in the second the velocity in one beam 
is different. The error velocity in the second case will be on average larger than 
the error velocity in the first case. The difference in the error velocities can bo 
from different currents or from errors caused by malfunctioning equipment. 

D. DEPLOYMENT AREA 

The RD Instuments ADCP remotely measures water flow velocity along the 
lines of position defined by the four narrow vertically inclined acoustic beams. 
To insure accurate current measurement it is necessary that the water mass in the 
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region of measuremenl be free of strong acoustic reflecting objects (e.g platform 
members, large cables, surface, etc.) within a +15 degrees conical sector along the 
direction of each of the beams. However, since only three beams are required for 
computation of three axes current components, in applications where potential 
interference objects may be close to one beam, the current vector may be calcu- 
lated from the other three beams. 

In shallow water vertical profiling applications, acoustic interference from the 
surface (upward looking ADCP) or bottom (downward looking ADCP) limit the 
vertical current measurement region to a maximum range defined by 

R{ max) = D X cos (j) (9) 



where 

• R(ma.x) = maximum profiling range, 

• D = distance to surface/bottom boundary, and 

• (f) = acoustic beam angle relative to the vertical. 

E. BOTTOM TRACKING OPERATION 

The ADCP can be used in a bottom tracking mode to give direct estimates 
of the velocity components of the ship. In the experimental data collected for this 
thesis we used the bottom tracking on for the runs 3 - 2, 2 - 3, and at the first part 
of the run 3-4. In order for the ADCP to gather bottom tracking data, bottom 
echoes must be distinguished from other echoes. The ADCP transmits a dedi- 
cated bottom track ping between current profiling pings. The number of profiling 
pings can be selected by the user. Consequently, the bottom tracking ping data 
is separated from the current profiling data. Bottom echoes are identified by vir- 
tue of their greater echo strength (the bottom returns a stronger echo than the 
echo from a profiling bin). 

Ship \elocity (relatit'e to the bottom) is measured in much the same way as 
current velocity (the Doppler shift of the backscattered bottom echo being pro- 
portional to the ship velocity). Bottom depth is detected by comparing recei\ed 
echo amplitude to a detection threshold which decreases with increasing range to 
the bottom. 
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F. VELOCITY PROFILES 



The most important feature of ADCP's is their ability to measure current 
profiles. ADCP's break up the velocity profile into uniform segments called depth 
cells. 

Each depth cell is similar to a single current meter, but there are two basic 
differences. The first difference is that the depth cells in an ADCP are alwa 3 's 
uniformly spaced while current meters can be spaced irregularly, the second is 
that the ADCP measures average velocity over the depth range of each depth cell 
while the current meter measures current only at the current meter Figure 7. 






y. Current Velocity Vector 

/I 

"<L:.^=5ic 

Avfio'jcj VerfQr^ W'lhin 
Complete Depth Cell Aieo 

; 

^ ^ 

! <! - 
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Figure 7. ADCP AND CURRENT ME! ERS 



The advantage of the ADCP averaging the velocity over the full range of a 
depth cell is that it reduces the effects of spatial aliasing (signals at frequencies 
higher than the time series can resolve arc mistaken for low frequency signals). 

G. ACOUSTIC DOPPLER CURRENT PROFILER DATA PROCESSING 
The raw data is collected in 128 four meter bins over a three minute ensemble 
interval. Generally speaking, either a MINI RANGER, LORAN C (calculated 
<j) , ?. from time differences or </> , 2 recorded from the display) or a GPS rccci\ cr 
are used for navigation and a position recorded every 30 seconds. By interpo- 
lation a position can be calculated at the end of each ensemble. 



The first steps in processing are the calculation of ship's velocity from the 
navigation data. From this navigation data the U and V components of ship's 
velocity are calculated. 

The next step in processing is the initial determination of the depth to which 
the data of each ensemble remains reliable. For our case the 7 , 8 , 9 bins were 
chosen. The basic criterion comes from the good percent of return echoes (BIN 
STATISTICS FILE). By subtracting the ship's velocity from the average velocity 
within the chosen reference layer an absolute reference layer velocity for each 
ensemble is obtained. The series of absolute reference velocities is then filtered 
with a low pass Hamming window filter. 

Once the absolute reference \ elocity is determined the velocity profiles of each 
ensemble with respect to the reference velocity are also determined, thus yielding 
the final profiles of absolute water velocity. The remaining profiles of absolute 
velocity are then av eraged over the time interval for each run. 

In appendix E there are 18 average current profiles from different navigation 
data. 



12 



IV. THEORY OF LORAN C 



A. INTRODUCTION 

LORAN C (LONG RANGE NAVIGATION) is a pulsed, low frequency 
(100 KHZ carrier) long-range hyperbolic navigation system. It operates on the 
principle that the difference in time of arrival of signals from two stations, ob- 
served at a point in the coverage area, is a measure of the difference in distance 
from the point of observation to each of the stations. 

Measurements of the ship's latitude and longitude ( (^ , / ) and Time Differ- 
ences (TD) from a LC 408 LORAN C receiver were recorded throughout the 
experiment (22 September 1989) at intervals of 30 seconds. 

B. SYSTEM DESCRIPTION 

LORAN C normally rcv-jUires signals from at least a Master and two Sec- 
ondary stations to gi\e a positional fix. 

Transmitters are grouped to form a "chain" of which one station is labelled 
the Master (M) and the others are called secondaries (X,Y,\V) as shown in Figure 

8 . 

The chain designators for LORAN C are 4 digit numbers which indicate the 
pulse Group Repetition Interval (GRI) in tens of micro-seconds. For example, the 
west coast LORAN C chain (Figure 8) is designated 9940 and has a GRI 99400 
micro-seconds ( see Figure 9) [LC 408, Operation Manual] 

C. PULSES-PHASE AND CODES-CYCLE SELECTION 

Each station of a LORAN C chain transmits groups of pulses (Table 3). 



Table 3. LORAN C PHASE CODES 



Master 


A + +-+-+- + 


B -h--h -h -f- -h -h - 


Secondary 


A -h-l- + -h-f---h 
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Figure 8. LOCATION OF WEST COAST LORAN C STATIONS 
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Figure 9. LORAN C GROUP REPETITION INTERVAL (GRI) AND TIME 
DIFFERENCES (TD) 



The Secondary transmits eight pulses and the Master nine (the last pulse fc>r 
indentification and blink alarm). The pulses are phase coded to improve the sig- 
nal to noise ratio through compression and to distinguish between Master and 
Secondary transmissions. 

Low frequency radio signals propagate over the earth's surface at nearly the 
velocity of light in a vacuum. However, Maxwell's equations dictate that the 
ground wave (surface wave) will be influenced slightly by the surface parametcis 
of geometry and electrical properties. 

In order to make the received 100 KHz signal more stable and reliable witliin 
a given coverage area, the LORAN C radio navigation system is designed as a 
pulse system which separates the ground wave from the sks^vave. 

Because the earth parameters remain nearly constant, LORAN C has dem- 
onstrated a repeatability of quite high accuracy. 

The high accuracy of LORAN C, despite long ranges from transmitters, is 
due to a technique called "cycle matching" [Bowditch, 1984]. 

The LC 408 tracks the third cycle cross-over path point which is very con- 
sistent between transmissions and less susceptible to skvxvave interference than 
later cycles. By tracking this cycle on all pulses, high accuracy is attained. Figure 
10 shows the third cycle tracking point with an example of sky wave interference. 
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D. TIME DIFFERENCE EQUATIONS 

The elapsed time between the arri\al of the signals from the Master and a 
Seeondary stations is called the "Time Difference" (TD). Each observed time 
difference, or rate, provides one hyperbolic line of position. By observing the 
transmissions from four stations, three hyperbolic rates are measured and the 
position ean be determined by either graphical or analytical techniques. 

The time differenee observed at a receiver is the difference in arrival times for 
signals from the master and one seeondary transmitter in the ehain. Because all 
transmitters share the same frequencies, their signals must be separated in time 
to prevent interference. 

The chain is syehronized so that the Master transmits first followed by each 
of the Secondaries. The transmission of each seeondary is specified by the emis- 
sion delay so that in the eoverage area signals from one station will overlap an- 
other [Schenebclc, 1979]. 

Suppose that at a point P there is a receiver. The observed time difference 
using the Master M and Secondaries X,Y,W are 
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Figure II. TYPICAL HYPERBOLIC FIX GEOMETRY 



TDW = ED^y+ tw- tsi (10) 

TDX = EDx+tx-t^j ( 11 ) 

TDY = EDy+ ty- (12) 
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where EDx, EDy are the emission delay for W, X, Y. The terms , ty 
are travel times from the W , X , Y to P. The term is the travel time from the 
Master M to P. 

In order to express the time difference as a function of geographic position , 
the travel time t is separated into additive terms 

r = -^xZ) + F (13) 



where: 

• C = free space propagation velocity, 

• n = index for refraction for a standard atmosphere, 

• D = geodetic distance from the transmitter to receiver, and 

• F = phase factor which corrects for effects of the earth's surface along the 
path. 

Substituting the last formula into the equations (10), (11), (12). gives the 
equations; 



rZ)l|-= £Z)„. + -^x (D,,,- ( 14 ) 

TD.\ = £D V + -^ X (£>,v - £.u) + Y - f 1/ ( 1 5) 

ro )■= £D X (£,.-£ „) + Fy - £,, ( 1 6) 

All these equations relate the time differences TDW.TDX.TDY to the dis- 
tances from each of the four transmitters. 

Since the latitude and longitude is mathematically calculated from the time 
difference numbers, the crossing angles, gradients, and signal to noise ratios are 
still important. 

E. HYPERBOLIC GRADIENT, CROSSING ANGLES 

For a gi\ en pair Master, Slave a family of hyperbolas separated by a costant 
\ alue can be plotted. 



IS 




The hyperbolas in Figures 12 and 13 are all separated by 10 microseconds 
[Gann]. That means that the distance on the baseline represents about 1852 
meters. The distance between successive hyperbolas increases as one moves to- 
wards the baseline extensions or away from the baseline. This change in the ac- 
curacy of a hyperbolic Line of Position (LOP) that occurs relative to the position 
of the LORAN C stations is known as the gradient. 

The angle between the LOPs is called the crossing angle Figure 13. The 
crossing angle must be bigger than 30 degrees and smaller than 150 degrees. This 
fact is common with any ranging positioning system, but is further limited for the 
LORAN C because the reception of three stations is required to define two LOPs 
(direct ranging navigation systems provide one LOP per station). 



19 




F. ERROR ANALYSIS FOR LORAN C 

An\" discussion of accuracy should begin with the definition of the types of 
accuracy. 

Repeatable accuracy (Repeatability, Precision) is the accuracy with which a 
navigator can return to a location, the coordinates of which have been previously 
measured using the same system. 

Random error is one which results from basic limitations in the method. The 
characteristics of this error can be determined by statistical analysis of a sufficient 
number of measurements. This type of error affects the repeatability. These er- 
rors are identified by their Gaussian distribution. 

Absolute accuracy is the accuracy with which a navigator can determine his 
point position in terms of a well defined coordinate system e.g,(</> , ). 
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Systematic errors result from a basic (but unrealized) fault in the method 
and cause the values to be consistently biased from the true value. It affects the 
absolute accuracy, and often cannot be detected by statistical analysis. 

The LC 408 microprocessor uses an exacting spherical trigonometric formula 
to derive its latitude and longitude positions. Also it is important to note that the 
LC 408 used for the calculations assumes that the LORAN signals travel over 
ideal all seawater paths from the transmitter to the receiver. Tp correct for path 
over land, the LC 408 allows use of "additional secondary factors"; this is typi- 
cally used only in small regions adjacent to known fixed locations. 

For the error analysis given shortly, the shape of the earth is approximated 
by a sphere, and it is assumed that all the west longitudes are positive . 

In order to find the standard deviations in the velocities, we established a 
fixed point in the roof of Spanagel Hall to enable an assessement of the standard 
deviations in the three time differences and after it to propagate the error at the 
positions and the velocities. 



Consider the spherical triangle PSM (see Figure 14 [Cross, 1981]). 




Figure 14. SPHERICAL TRIANGLE SPANAGEL-M-NORTH POLE 
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Where (M) can be the Master, S can be the point on the Spanagel roof, and 
P can be the North pole. We try to find an expression for the distance between 
Master and Spanagel roof with respect lo <{) , A (coordinates of the point on the 
Spanagel rooQ, 



m = 90 — (f) (17) 

s = 90-4>m ( 18 ) 

( 19 ) 

where 4>xf and are the coordinates of the Master station. Using the law of 
cosines for the side p, 

cos p = cos m cos s + sin m sin s cos P (20) 

Substitute (17) , (18) , (19) into (20) to yield 

cos p = sin sin 0 + cos cos ({) cos(2 — /^.^^/) (2 1 ) 

Differentiate each part of (21) 

— sin pdp = cos (/) V/ sin 4>d(f) \^ + sin (/> v/ cos 4)d(j) 

— sin y cos 4> cos{/. — }.^j)d4)^{ 

— cos 0,\/sin (f) cos(/ — A^{)d4) (22) 

— cos 0 cos 4> sin(/ — /.\t)d?. 

+ cos cos (f) sin(/ — A^,i)d/.^j 

Since the position of Master station is a fixed point d/.^ = <70,v/ = 0 equation (22) 
becomes 



- sin pdp = [ sin 4>m cos 4> — cos sin 4> cos(/ - 

— cos <p ^ cos 4> sin(/ — ^. \i)d/. 



(23) 



Wc now introduce into each term in the right part of (23) the factor sin p. Using 
the relationship between two angles and three sides 
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sin p cos S = cos s sin m — sin s cos tn cos P (24) 

Substituting for s, m and P as appropriate, we have 

sin </) cos 4> — cos 4> :^^ sin 4> cos(/ - = 

cos(90 - sin(90 - <j>) — sin(90 - ^ v/) cos(90 - 4>) cos(2 - ?.^) = (25) 

sin p cos S 

From the law of sines we know that 



sin m _ sin p 
sin M sin P 



(26) 



OR 

sin m sin P = sin /? sin M (2S) 

Substituting for m and P and multiplying both sides by cos , we have 
cos 0 )/ cos (j) sin(/ - 2 j/) = 

cos sin(90 — 4>) sin(2 — = (29) 

sin p sin M cos (/> ^^ 

Substituting (25) and (29) into (23). we have 

dp = — cos Sd(p + sin M cos 4> \idA (30) 

Actually in order to cont ert the angle dp into distance da (arc length), using 
spherical approximation multiply by the mean radius of the earth R 

da = — R cos Sd(f) + R sin M cos <p^\fd/. (3 1 ) 

From the last equation we have 

4^=-RcosS (32) 

C(p 
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= R cos <f)>^ sin M 



( 33 ) 



where M and S are the interior angles of the triangle PSM. In similar way 
form the equations for the other spherical triangles 



da^ 



— R cos Sj 



dai 

~dl 



= — Rcos (/jji'Sin X 



The sign changes in the second equation due to different geometry 



da-, 

= — ^ cos S2 



cai 

= R cos </)ysin )' 



ca^ 

ttt = - Rcos 5, 

C(p 



can 

= R cos (hu-sin ir 
c/. 



Next let 



= a{(j), /.) - aj(0, ;.) 



di = a{(}), ;.) - a 2 {(f), /) 



= a{(f), ;.) - a^icf), /) 

Suppose that the velocity of the the electromagnetic propagation is constant 
equal (C) then: 



can 

( 34 ) 

( 35 ) 

( 36 ) 

( 37 ) 

( 38 ) 

( 39 ) 

( 40 ) 

( 41 ) 

( 42 ) 
and 
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(43) 



A 



'1 



A 

C 



A 



'2 




(44) 




(45) 



We propagate the errors into these expressions by forming the Jacobian ma- 
trix. We assume here that time is subject to negligible error. This matrix is given 
by 




c’<f] 

cd. 

C(p 

ddi 



dd] 

c?. 

dd. 




(46) 



and 

^d^ ca{<j),?d ca^{4),/.) 

d(j) c4> 

= — R cos S + Rcos Si (47) 

= R{ cos Sj — cos 5) 

cd2 da{(j), /) ca2(4>, /) 

d<f) d(j) d(f) 

= — R cos S + R cos S 2 (48) 

= R{ cos S 2 — cos S) 

ddj ca{(j). /) da 2 {<l>- 

C(p C<j) C(j) 

= — R cos S + R cos S 3 (49) 

= R[ cos S 3 — cos S) 
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( 50 ) 



dd^ da{(^,X) dai{<f),).) 

dA d/. 5/ 

= R cos 4>f^ sin M + R cos 4>;^ sin X 

= R{ cos (^jV/ sin M + cos (f>x sin X) 

ddi da{(f), /) da2i<j>, A) 

Sa dA dA. 

= R cos 0 ;v/ sin AI — R cos 0 y sin K (51) 

= R{ cos sin M — cos ^y^sin F) 

ddj da{(f>, /) da^i^), a) 

dA. 6 a. dA. 

= R cos sin M — R cos (f)iySin IV (52) 

= /?( cos y/ sin M — cos <piy sin 11") 

C = 3 X 10^('\// sec) = 3 X \0~\Mlnsec) 



Now define the variance covariance matrix of LORAN time differences, as- 
suming that there is no correlation between the three time differences. 



V _ 




0 

0 



0 

0 



1 




(54) 



The LC 408 is provided with control of the filtering time constant. Selectable 
time constants are 0 seconds for no filtering, 5 seconds for good response, 10 
seconds for slower response, and 20 seconds for slow craft or monitor applica- 
tions. 

In order to measure the standard de\iations for the time differences, we re- 
corded data for four days on Spanagel Roof, each day (24 hours) with a different 
time constant. This resulted in four different variance covariance matrices. 
There are all gi\ en in Appendix A. 

From these results we conclude that the filters don't work for the time dif- 
ference between Master and Secondary Whisky because the standard deviation 
remains approximately the same for the different time filters. We ha^e no ex- 
planation for this beyond it being an instrumental problem. We also conclude 
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that for monitoring applications, the longer the filtering time, the smaller the 
standard deviations. 

The weight matrix P — L"' is given by 



P = 




0 



1/(7^ 

0 



0 

0 




(55) 



The variance covariance matrix for the position [Uotila, 1986] will be: 

1(0, /) = (y^x P X (56) 



Define the matrix Zq as 



V 

-o - 



'0 

0 




(57) 



Now we form the block matrix Z ( 0, /, which is the t ariance 

covariance matrix for two positions (the error characteristcs are the same for the 
two positions). If it is assumed that there is no correlation between the two posi- 
tions, then 



— (01 7 -], 02 /- 2 ) 



!(</>!, 7.,) 
^0 



V 

^0 






(58) 



In order to con\ert the difference in latitude and difference in longitude (radians) 
between two positions in meters, apply 

A = A<f)x R (59) 

B = Aax cos (^^x R (60) 

where 0,„ is the mean latitude of the work area, and R is the mean radius of the 
earth R = 6371 KM. Let 



C, = 7? 



(61) 
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C2 = COS X 



(62) 



The components of the ship velocity V , U (Northing, Easting) are then given by: 





(02 - 0l) 



(63) 






At 



£i 

At 



X i^-2 - /l) 



The Jacobian matrix for the velocity will be 



./ = 




dV 

dh 

du 

dh 



dV dV 

O/j C02 0/2 

dU dU eu 

d/.y 84>2 C/-2 



and 




-Cl 

0 



0 Cl 0 ■ 

“C2 0 C2 



(64) 



(65) 



( 66 ) 



Consider time inter\’als between positions 3 , 6 , 12 , 20 and 40 minutes. The 
\ ariance co\ ariance matrix for the \elocity components will be: 



V 

u 



= 7x1 



(?), 2 |. 4 > 2 . /-2 



X J‘ 



(67) 



Appendix B gives 5 tables for standard deviations and covariances of t elocity 
components with different time intervals between positions. From these tables 
we see clearly that as the time interval between the positions increases the stand- 
ard deviations for the V and U components decreases (assuming constant course 
and speed of the ship for the gi\en interval). Also we note that when time aver- 
aging over 20 minutes or longer, the choice of filter is not critical. 

Appendix C gi\'es a program (DRIVER FORTRAN) written by the NOAA 
which calculates the geographic position of the ship from two LORAN Time 
Differences (TD) and a Dead Reckoning (DR) position of the ship (iinerse com- 
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putation). It can also be used to convert the geographic position of the ship into 
LORAN Time Differences (forward computation). 
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V. THEORY OF MINI RANGER 



A. INTRODUCTION 

The Mini Ranger Falcon was used during the cruise for navigational posi- 
tioning. The standard Mini Ranger operates at up to 37 kilometers (about 20 
nautical miles) with a probable range measurement error of 2 meters. 

It operates at microwave frequencies and requires that line of sight be main- 
tained between the reference stations and the receiver transmitter. Significant 
obstuctions such as land masses, buildings, or dense foliage will interfere the op- 
eration of the system. 

B. REFERENCE STATIONS 

The positions of the reference station sites in the UTM coordinate system are 
[Krioneritis. 1989] listed in Table 4. 



Table 4. REFERENCE STATION SITES 



STATION 


X EASTING 


Y NORTHING 


CODE 


TREVOR 


585260.161 M 


4092490.284 M 


15 


PACK 


609863.128 M 


4076611.345 M 


1 


HAYES 


607621.289 M 


4055915.264 M 


12 



C. RANGE POSITIONING 

The Mini Ranger determines a two range position when the lengths of all 
three sides of a triangle are known (trilateration). Suppose that a certain task re- 
quires acti\ity between points C , D , E (Figure 15), the reference stations arc at 
points A and B. in the acceptable area the angle of intersection between the two 
range lines is between 30 and 150 degrees. 
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Figure 15. PLANNING THE MOST DESIRABLE SITE 



Arcs AXB and AYB define the maximum and minimum distance to maintain 
between 30 and 150 degrees respectively. 

Generally the Mini Ranger uses the method of least squares to calculate po- 
sitions from three or more reference stations. Operation in areas where the gc- 
t ometry is poor will degrade positional accuracy [Mini Ranger operation manual, 
1981]. 

D. ERROR ANALYSIS FOR MINI RANGER 

While the Mini Ranger measures ranges well, the actual position error re- 
sulting from a set of range readings depends upon a number of contributing fac- 
tors: 

1 . system error, and 

2. geometry changes. 
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Figure 16 shows the maximum positional error for two range geometries of 30 , 
90 and 150 degrees [Mini Ranger operation manual, 1981]. 



POSITION ERROR = 



STANDARD ERROR 
sin 6 



(68) 



The position error for various two range geometries, assuming range error of a 
single measurement of 2 meters are given in Table 5. 

Table 5. POSITION ERROR FOR DIFFERENT 



CROSS ANGLES 



CROSSING ANGLE 


POSITION ERROR 


150 degrees 


7.7 meters 


120 degrees 


4.0 meters 


90 degrees 


2.8 meters 


60 degrees 


4.0 meters 


30 degrees 


7.7 meters 



For this experiment we'll take the information for the variance covariance 
matrix of the positions from the Mini Ranger data processing (see, for example, 
Krioneritis, 1989). Define the matrix Zo as 



We now form the block matrix L ( which is the \'ariance 

coN’ariance matrix for two positions. If it is assumed there is no correlation be- 
tween the two positions, then 




(69) 



The variance co\'ariance matrix for the position will be 





(70) 
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• t 2 meteis 





Figure 16. MAXIMUM POSITION ERROR 



uxo'i/iyi) 



I(Ai, >'i) 

V 

^0 



2^0 

2(a’ 2, 1 y 



Define the velocity components as 



V = 



y2 - > 1 

At 



(71) 



(72) 



U = 




At 



(73) 
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As we did with LORAN C, we propagate errors into these expressions by forming 
the Jacobian matrix for the velocity components. We again assume that time is 
subject to negligible error. This matrix is given by 




IK 

dX\ 

ILL 

dXi 



8V dV 

8X2 

8U 8U 

ai'i 8X2 



IK 

8Y2 

ILL 

8Y2 



OR 





-1 0 n 

0 1 oj 



The variance covariance matrix for the velocity components will be 



(74) 



(76) 



(77) 



Table 6. MAX-MIN VALUES OF STANDARD DEVIATIONS AND 
COVARIANCES FOR VELOCITY COMPONENTS FOR 3 MINUTE 
TIME INTERVALS BETWEEN POSITIONS (MINI RANGER) 



MAX-MIN 

VALUES 


Oy 




^ V, V 


MAXIMUM 


1.6 cm; sec 


1.8 cm, sec 


l.ScnY/ sec- 


MINIMUM 


1 .2 cm sec 


1 .3 cm sec 


Q.lcnYj scc‘ 



Table (6) shows that the the standard error in the velocity components 
doesn't \ ary greatly gi\en the geometry in the work area. The results come from 
6 different points in the whole work area. 

Appendix D gi\es the \ elocity error results for different time intervals. Simi- 
lar to the LORAN (Appendix B), the standard de\iations for the velocity com- 
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ponents decrease when the time interval between the positions increases 
(assuming constant course and speed of the ship for the given time interval). 

Comparing the theoretical results from Appendices B and D, we can see that 
the MINI RANGER is not only a more stable navigation system than the 
LORAN C, but also that it provides in 3 minutes the level of accuracy in the U 
and V velocity components which require 20 minutes to obtain with LORAN C. 
Because of the practical difficulties in maintaining a constant ship's heading and 
speed for 20 minutes, and thus the need for an averaging process over such a time 
interval, the value of using a MINI RANGER wherever possible is even further 
enhanced. 
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VI. RESULTS 



In this chapter there are eleven Tables which summarise the results of the 
computation of the ship velocity components from MINI RANGER, LORAN 
Time Differences (TD), LORAN (DISPLAY) and Bottom Tracking. We choose 
to separate the LORAN (TD), and LORAN (DISPLAY) results because of the 
rounding off which occurs w’hen using the LORAN (DISPLAY). 

The LORAN results given here were obtained using different filters (2 - I, 1 
- 3 filter with time constant 5 seconds, 1 - 2 filter with time constant 20 seconds, 
3 - 2, 2 - 3 filter with time constant 10 seconds, 3-4 filter with time constant 0 
seconds ). 



Table 7. SHIP VELOCITY COMPONENTS FOR MINI RANGER, LORAN 
(TD), LORAN (DISPLAY! FROM POINT I TO POINT 2 AT 3 MIN- 
UTE TIME INTERVAL. 



TIME 


MINI RANGER 


LORAN (TD) 


LORAN (DIS- 
PLAY) 


U 

(cm sec) 


V 

(cm sec) 


U 

(cm sec) 


V 

(cm sec) 


U 

(cm sec) 


V 

(cm sec) 


14 24 1 1 


-126.5 


-277.1 


-125.4 


-273.5 


-126.2 


-263.3 


14 27 10 


-125.9 


-276.0 


-123.2 


-274.8 


-124.0 


-286.4 


14 30 10 


-123.3 


-277.6 


-120.6 


-278.1 


-115.7 


-271.5 


14 33 1 1 


-124.4 


-280.5 


-124.4 


-274.6 


-124.0 


-282.8 


14 36 11 


-106.3 


-282.8 


-107.4 


-282.5 


-109.3 


-281.6 


14 39 1 1 


-108.7 


-285.1 


-103.0 


-296.2 


-105.2 


-275.6 


14 42 10 


-109.3 


-280.7 


-116.5 


-284.5 


-113.4 


-268.4 
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Table 8. SHIP VELOCITY COMPONENTS FOR MINI RANGER, LORAN 
(TD), LORAN (DISPLAY) FROM POINT 2 TO POINT 1 AT 3 MIN- 
UTE TIME INTERVAL 



TIME 


MINI RANGER 


LORAN (TD) 


LORAN (DIS- 
PLAY) 


U 

(cm sec) 


V 

(cm/Sec) 


U 

(cm/sec) 


V 

(cm/sec) 


U 

(cm ''sec) 


V 

(cm, sec) 


14 57 10 


135.0 


265.3 


134.3 


254.5 


135.0 


258.0 


15 00 10 


133.7 


261.3 


134.4 


256.1 


132.2 


240.7 


15 03 10 


13S.9 


255.4 


134.7 


253.7 


138.1 


265.0 


15 06 10 


154.4 


252.0 


154.0 


245.6 


148.7 


244.1 


15 09 10 


155.9 


253.1 


154.7 


256.5 


151.7 


247.9 


15 12 10 


148.8 


255.4 


146.2 


252.3 


154.1 


261.3 


15 15 10 


155.7 


251.8 


154.9 


250.7 


143.2 


237.2 



From Tables 7 and S we can sec a discontinuity beginning at 14 36 11 near 
the end of leg 1 - 2 and ending at 15 03 10 at the beginning of leg 2 - 1. This 
discontinuity occurs at the same position down each leg and is of approximately 
the same magnitude in each case. While it is tempting to view this discontinuity 
as being due to a surface current, the inconsistency in sign prohibits such an ex- 
planation. It can only be concluded therefore that it is due to a ship navigation 
correction. 

In these tables we see clearly that the LORAN (DISPLAY) results show a 
much gcater variability than the LORAN (TD). Also if we calculate the mean 
value and the standard deviation of the mean of the absolute difference of the U 
and V components between the MINI RANGER and LORAN (TD), and the 
MINI RANGER and LORAN (DISPLAY), we have Table 9. 
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Table 9. COMPARISON BETWEEN MINI RANGER AND LORAN (TD) AND 



THE MINI RANGER AND LORAN (DISPLAY) FOR LEGS 1 - 2 AND 
2 - 1. 



LEG 


MINI RANGER - LORAN (JD) 


MINI RANGER -LORAN (DISPLAY) 




A V (a^y) 


A U (cr^j,0 


A V (<T^v) 


1 - 2 


2.9(0,9) 


4.5(2. 1) 


2.9(0.9) 


7.9(1.8) 


2 - 1 


1. 5(0.5) 


4.5(1.2) 


4.2(1.6) 


10.2(2.1) 



After using the FISHER-BEHERENS test [Hamilton, 1964] to compare AU from 
MINI RANGER-LORAN (TD) against AU from MINI RANGER-LORAN 
(DISPLAY) and AV against AV we see that these are statistically different at 
95% confidence interval in 3 out of 4 cases. This provides strong evidence to 
suggest that the results from LORAN (TD) and LORAN (DISPLAY) are from 
different statstical populations. 



Table 10. SHIP \ ELOCITY COMPONENTS FOR MINI RANGER, LORAN 
(TD), LORAN (DISPLAY’) FROM POINT 1 TO POINT 3 AT 3 
MINUTE TIME INTERVAL 



TIME 


MINI RANGER 


LORAN (TD) 


LORAN (DIS- 
PLAY) 


U 

(cm sec) 


V 

(cm sec) 


U 

(cm/sec) 


V 

(cni;Sec) 


U 

(cm sec) 


V 

(cm sec) 


15 30 10 


-272.4 


-116.1 


-272.0 


-118.6 


-272.8 


-123.8 


15 33 1 1 


-273.1 


-115.8 


-274.6 


-116.5 


-253.5 


-106.8 


15 36 11 


-277.4 


-117.8 


-276.6 


-116.1 


-291.8 


-126.8 


15 39 10 


-277.1 


-117.9 


-276.8 


-121.4 


-273.7 


-104.1 


15 42 10 


-268.4 


-122.7 


-269.5 


-123.2 


-274.8 


-130.6 


15 45 10 


-266.7 


-126.4 


-268.6 


-130.2 


-264.4 


-134.0 
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15 48 10 


-265.9 


-123.9 


-264.7 


-119.4 


-264.7 


-123.9 


15 51 10 


-268.9 


-126.7 


-265.9 


-119.8 


-264.3 


-113.5 


15 54 10 


-269.8 


-123.9 


-263.8 


-118.7 


-276.3 


-123.9 


15 57 10 


-266.7 


-126.1 


-267.2 


-124.8 


-250.8 


-125.8 


16 00 10 


-269.5 


-120.5 


-269.5 


-114.8 


-269.9 


-101.3 


16 03 11 


-265.9 


-127.2 


-259.5 


-127.9 


-273.8 


-140.1 


16 06 10 


-268.6 


-125.5 


-267.4 


-125.2 


-260.6 


-123.8 



Table 11. COMPARISON BETWEEN MINI RANGER AND LORAN (TD) AND 
THE MINI RANGER AND LORAN (DISPLAY') FOR LEG 1 - 3. 



LEG 


.MINI FLANGER - LOR.AN' (YD) 


.MINI R.ANGER -LOR.AN (DISPLAY) 


V (<7jt ) 


A V (c-^, ) 




A \ 


1 - 3 


L9(0.6) 


2.9(0.6) 


7.0(1.7J 


7.9(1.-) 



Tabic 10 shows similar results for leg 1 - 3. From Table 1 1 it is also obt ious 
that the results which comes from the LORAN (TD) and the LORAN (DIS- 
PLAY') arc different for leg 1 - 3. 



Table 12. SHIP VELOCITY COMPONENTS FOR MINI RANGER. LORAN 
(TD), LORAN (DISPLAY^, BOTTOM TRACKING FROM POINT 3 TO 
POINT 2 AT 3 MINUTE TIME INTERVAL 



TIME 


MINI R.ANGER 


LOR.\N* 0 D) 


LOILAN (DIS- 
PLAY) 


BOETO.M 

TR.ACKING 


L 

(cm sec) 


(cm sec) 


U 

(cm sec) 


\’ 

(cm sec) 


U 

(cm sec) 


(cm sec) 


U 

(cm sec) 


\* 

(cm sec) 


16 16 18 


272.4 


-103.5 


267.9 


-104.4 


259.7 


-105.4 


263.5 


-116.5 


16 19 18 


277.9 


-93.2 


280.2 


-94.6 


274.1 


-84.0 


2"1.4 


-103.5 


16 22 IS 


275.5 


-95.3 


272.1 


-95.5 


279.7 


-111.8 


2()S.5 


-106.0 


16 25 18 


277.3 


-99.2 


275.1 


-97.5 


276.6 


-80.7 


2“'3.3 


-110. 1 
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16 28 17 


278.0 


-98.2 


282.2 


-95.4 


275.7 


-98.7 


272.4 


-108.1 


16 31 19 


280.1 


-106.6 


279.8 


-97.7 


281.3 


-99.1 


274.2 


-109.5 



Table 13. SHIP VELOCITY COMPONENTS FOR MINI RANGER, LORAN 
(TD), LORAN (DISPLA’VO, BOTTOM TRACKING FROM POINT 2 TO 
POINT 3 AT 3 MINUTE TIME INTERVAL 



TIME 


MINI RANGER 


LORAN (TD) 


LORAN (DIS- 
PLAY) 


BOTTOM 

TRACKING 


L 

(cm sec) 


(cm sec) 


U 

(cm sec) 


V 

(cm sec) 


U 

(cm sec) 


V 

(cm sec) 


L 

(cm sec) 


V 

(cm sec) 


16 52 18 


-260.3 


124.3 


-256.1 


124.3 


-256.5 


124.3 


-253.0 


128.5 


16 55 17 


-259.2 


118.4 


-261.1 


114.4 


-271.0 


121.5 


-252.8 


124.8 


16 58 17 


-262.8 


116.4 


-264.7 


115.9 


-249.9 


108.1 


-257.1 


125.9 


17 01 17 


-269.5 


122.7 


-266.1 


118.0 


-277.1 


126.5 


-262.8 


128.4 


17 04 18 


-271.1 


117.8 


-268.9 


120.9 


-266.6 


118.5 


-264.9 


125.1 


17 07 18 


-264.2 


116.4 


-265.3 


108.3 


-266.8 


114.2 


-258.8 


123.0 


1" 10 IS 


-267.0 


115.4 


-263.6 


119.4 


-260.2 


120.8 


-262.9 


123.2 


17 13 17 


-264.0 


121.9 


-260.9 


120.7 


-264.3 


116.0 


-258.4 


128.0 



From Tables 12 and 13 we can see that the more consistent results comes 
from MINI RANGER. The variability of the LORAN (TD) results are smaller 
than the variability of the LORAN (DISPLAY) results. We also see that the 
Bottom Tracking, while gi\ing values to the same internal consistency as the 
MINI RANGER, gives results which are systematically biased. 



Table 14. COMPARISON BETWEEN MINI RANGER AND LORAN (TD) AND 
THE MINI RANGER AND LORAN (DISPLAY^ FOR LEGS 3 - 2 AND 



2 - 3. 



LEG 


MINI ILANGER - LOILAN 
HD) 


MINI RANGER -LOR.AN 
(DISPLAY) 


.MINI R.ANGER ■ BOT- 
TOM IR.ACKING 


A U (o-^f ) 


A V (a^y) 


A U 


A V 


A U ((7,,.) 


A \- 
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3 - 2 


2.8(0.6) 


2.6(1.3) 


4.1(1.8) 


9.1(3.I) 


6.3(0.7) 


9.6(1. 4) 


2 - 3 


2.6(0.4) 


3.2(0.9) 


6.3(1.6) 


3.7(1.0) 


5.9(0.4) 


6.7(0.6) 



From Table 14, after using the FISHER-BEHERENS test to compare AU 
from MINI RANGER-LORAN (TD) against AU from MINI 
RANGER-LORAN (DISPLAY) and AV against AV we see that these are sta- 
tistically different at a 95% confidence interval in 2 out of 4 cases. 



Table 15. SHIP VELOCITY COMPONENTS FOR MINI RANGER, LORAN 
(TD), LORAN (DISPLAVl. BOTTOM TRACKING FROM POINT 3 TO 
POINT 4 AT 3 MINUTE TIME INTERVAL 



TI.ME 


MINI R.'WGER 


LOR.AN (TD) 


LORAN (DIS- 
PLAY) 


BOTTOM 

TFLACKING 


L 

(cm sec) 


\’ 

(cm sec) 


U 

(cm sec) 


(cm sec) 


U 

(cm sec) 


(cm sec) 


U 

(cm sec) 


(cm sec) 


17 22 17 


-9.4 


2S6.1 


-12.1 


283.2 


-10.2 


275.7 


0.5 


282.2 


17 25 17 


-11.7 


288.0 


-12.8 


284.2 


-10.6 


288.9 


-1.4 


285.9 


17 28 17 


-12.5 


289.3 


-9.8 


285.3 


-10.9 


287.5 


-4.1 


2S7.0 


17 31 17 


-4.9 


295.7 


-8.7 


295.9 


-15.1 


293.5 


1.5 


290.6 


17 34 17 


-8.7 


293.5 


-11.7 


292.8 


-7.1 


293.0 


-1.9 


288.4 


17 37 19 


-19.4 


293.5 


-18.3 


286.7 


-23.6 


274.5 


-11.8 


290.1 


17 40 IS 


-21.6 


297.3 


-20.0 


296.2 


-20.0 


307.0 


-15.1 


293.2 


17 43 IS 


-16.2 


288.6 


-15.9 


2SS.8 


-7.9 


288.8 


-9.2 


286.6 


17 46 18 


-24.9 


292.2 


-26.4 


290.5 


-29.5 


299.3 


-20.2 


287.9 


17 49 18 


-37.5 


290.6 


-42.0 


284.7 


-41.6 


278.8 


-33.7 


286.8 


17 52 18 


-32.2 


288.7 


-34.5 


291.0 


-41.3 


287.5 


-27.7 


286.4 


17 55 17 


-28.0 


286.0 


-29.9 


281.1 


-23.4 


287.0 


-21.4 


282.5 


17 58 18 


-22.2 


281.3 


-13.6 


286.1 


-16.6 


247.7 


-15.3 


2’8.6 


18 01 18 


-22.3 


281.0 


-23.4 


276.0 


-25.7 


2S0.0 


-18.5 


276.7 


18 04 18 


-25.3 


275.4 


-26.1 


275.1 


-20.0 


283.2 


-17.:* 


273.9 


IS 07 IS 


-23.4 


274.6 


-24.9 




-28.3 


266.3 


-18.2 


2"1.6 
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18 10 18 


- 22.3 


271.0 


- 19.2 


273.1 


- 22.3 


258.7 


- 14.8 


268.9 


18 13 18 


- 28.7 


267.2 


- 30.2 


267.5 


- 23.8 


285.7 


- 23.1 


266.0 


18 16 18 


- 29.8 


262.8 


- 37.4 


257.3 


- 37.7 


262.5 


- 21.4 


265.1 


18 19 18 


- 30.9 


264.3 


- 20.7 


269.9 


- 20.0 


268.3 


- 17.1 


262.9 


18 22 47 


- 25.8 


261.2 


- 33.0 


257.9 


- 32.3 


256.5 


- 


- 


18 25 47 


- 31.3 


256.1 


- 37.8 


253.8 


- 41.2 


247.6 


- 


- 


18 28 47 


- 26.8 


266.0 


- 29.1 


262.2 


- 26.1 


262.7 


- 


- 


18 31 47 


- 29.1 


263.0 


- 25.7 


264.4 


- 28.0 


251.7 


- 


- 


18 34 47 


- 29.1 


264.7 


- 30.2 


263.2 


- 24.5 


281.2 


- 


- 


18 37 47 


- 27.5 


263.9 


- 27.2 


262.0 


- 24.9 


245.9 


i , 

•tJ 


- 


18 40 47 


- 33.6 


261.1 


- 35.8 


257.6 


- 32.8 


264.4 


- 


- 


18 43 46 


- 39.4 


256.5 


- 38.2 


261.5 


- 32.5 


251.1 


- 


- 


18 46 46 


- 32.9 


259.8 


- 34.4 


256.8 


- 45.4 


269.8 


- 


- 


18 49 46 


- 41.9 


257.5 


- 40.8 


258.9 


- 41.9 


258.0 


- 


- 


18 52 46 


- 43.3 


256.0 


- 43.3 


257.9 


- 44.9 


262.8 


- 


- 


18 55 47 


- 37.7 


257.9 


- 36.2 


255.3 


- 33.2 


248.9 


- 


- 


18 58 47 


- 32.5 


255.4 


- 31.0 


259.6 


- 32.5 


262.2 


- 


- 


19 01 47 


- 33.9 


252.4 


- 32.8 


251.9 


- 20.7 


231.4 


- 


- 


19 04 47 


- 32.4 


255.7 


- 32.4 


253.1 


- 49.4 


274.3 


- 


- 


19 07 46 


- 32.1 


253.1 


- 32.1 


253.1 


- 28.7 


237.3 


- 


- 


19 10 46 


- 43.5 


253.2 


- 41.2 


255.3 


- 41.6 


273.4 


- 


- 


19 13 4 " 


- 4“.8 


256.5 


- 48.1 


253.2 


- 49.3 


252.9 


- 


- 


19 16 46 


- 53.7 


255.4 


- 56.4 


251.4 


- 41.2 


232.0 


- 


- 


19 19 46 


- 5 ".l 


256.5 


- 56.3 


253.0 


- 65.8 


273.3 


- 


- 


19 22 47 


- 46.2 


259.8 


- 51.9 


254.7 


- 45.5 


242.7 


- 


- 


19 25 46 


- 42.0 


263.6 


- 45.8 


259.4 


- 53.3 


257.7 


- 


- 


19 28 46 


- 4". 7 


268.0 


- 49.6 


264.9 


- 49.6 


274.2 


- 


- 


19 31 46 


- 50.9 


270.9 


- 57.6 


263.3 


- 49.7 


262.8 


- 


- 


19 34 47 


- 59.5 


267.2 


- 60.6 


262.5 


- 65.9 


258.0 


- 


- 


19 3 ' 46 


- 54.9 


268.8 


- 59.0 


262.9 


- 53.4 


252.9 


- 


- 
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19 40 47 


-52.7 


272.5 


-51.9 


268.0 


-65.9 


283.8 


- 


- 


19 43 47 


-60.4 


276.0 


-66.0 


270.3 


-66.0 


273.6 


- 


- 


19 46 46 


-57.9 


279.1 


-64.4 


272.7 


-61.7 


273.4 


- 


- 


19 49 46 


-45.2 


273.9 


-52.7 


268.4 


-52.3 


264.4 


- 


- 


19 52 47 


-39.1 


276.8 


-36.9 


266.9 


-34.2 


266.9 


- 


- 


19 55 47 


-21.1 


278.8 


-23.0 


275.8 


-28.6 


273.6 


- 


- 


19 58 46 


-21.6 


281.4 


-25.0 


276.6 


-16.3 


273.3 


- 


- 


20 01 47 


-21.8 


276.8 


-26.7 


274.2 


-33.1 


278.7 


- 


- 


20 04 47 


-19.5 


277.7 


-22.2 


271.3 


-20.7 


268.3 


- 


- 


20 07 47 


-26.8 


279.4 


-30.5 


269.9 


-24.5 


268.5 


- 


- 


20 10 46 


-23.4 


274.2 


-24.6 


274.9 


-25.0 


268.5 


- 


- 



Table 16. COMPARISON BETWEEN MINI RANGER AND LORAN (TD) AND 
THE MINI RANGER AND LORAN (DISPLAY) FOR LEG 3 - 4. 



LEG 


.MINI RANGER - LORAN 
TD) 


MINI RANGER -LORAN 
(DISPLAY) 


MINI IL^NGER - BOT- 
TOM TR.ACKING 




A V (-7^ ) 


A U 


A V (<7^^.) 


A G (-7^^,.) 


A V (<7^,,) 


3 - 4 


2.S(0.3) 


3.5(0.3) 


5. 1(0.5) 


8.6(0.8) 


7. 1(0.5) 


3.K0.3) 



Table 15 again yields similar results for leg 3 - 4. From Table 16 it is again 
obvious that the results which comes from the LORAN (TD) and the LORAN 
(DISPLAY) are different. 

In Table 17 are the average velocity components for the different positioning 
systems and for the bottom tracking for each leg. From this table we can see; 

• Data averaged over times of about 18 minutes (minimum) show no statis- 
tically significant difference between MINI RANGER. LORAN (TD) and 
LORAN (DISPLAY). Any one of these positioning systems can be used for 
these averaged times. This is due to the fact that the variations of the ship 
course and speed in the aN eraged period of time completely swamp the dif- 
ferent precision which the different navigation systems have in their calcu- 
lation of the U and V components of the ship. The minimum time interval 
from which the averaged data over the different navigation systems become 
the same statistically depends upon the ship's variations in course and speed. 
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• The results from Bottom Tracking are clearly systematically different from 
those of the others systems. They do, however, show the same level of 
internal consistency as both the MINI RANGER and the LORAN (TD) 
results. 

• The much larger standard deviations on the LORAN (DISPLAY) results 
are indicative of the much greater variability in computed U and V compo- 
nents when using this system. 



Table 17, AVERAGE VELOCITY COMPONENTS FOR MINI RANGER, 



LORAN (TD), LORAN (DISPLAY), BOTTOM TRACKING 





MINI RANGER 


LORAN (TD) 


LORAN (DIS- 
PLAY) 


BOTTOM TRACK- 
ING 


U (<7l) 


V (<7v) 


t.' (c-t,) 




f (cr,-) 


\- (ay) 


(cf) 


V (<7y) 


1 - 
2 


-117.7(3.4) 


-279.9(1.2) 


-117.2(3.3) 


-276.7(2.1) 


-116.8(3.1) 


-275.6(3.2) 


- 


- 


2 . 
1 


146.1(3.7) 


256.3(1.9) 


144.7(3.8) 


252.7(1.4) 


143.3(3.2) 


250.6(4.1) 


- 


- 


1 - 
3 


-270.1(1.1) 


-122.3(1.1) 


-268.9(1.4) 


-121.2(1.3) 


-268.5(2.9) 


-121.4(3.2) 


- 


- 


3 - 
2 


276.8(1.1) 


-99.3(2.1) 


276.2(2.2) 


-97.5(1.5) 


274.5(3.1) 


-96.6(4.9) 


270.5(1.6) 


-108.9(1.8: 


3 


-264.7(1.4) 


119.1(1.1) 


-263.3(1.4) 


117.7(1.7) 


-264.1(3.0) 


118.7(2.1) 


-258.8(1.6) 


125.8(0.8) 


3 . 

4 


-32.5(1.8) 


271.1(1.7) 


-33.8(1.9) 


269.1(1.7) 


-33.6(2.1) 


268.7(2.1) 


- 


- 


3 - 
4A 


-21.6(1.9) 


283.3(2.4) 


-21.8(2.1) 


281.8(2.3) 


-21.8(2.3) 


282.1(2.7) 


-14.5(2.1) 


280.6(2.1) 



In Appendix E there are 12 figures which show the average current profiles 
from point 1 - 3, 3 - 2, 2 - 3, 3 - 4 with the three different navigation data from 
MINI RANGER, LORAN (TD), LORAN (DISPLAY). These plots show us 
the same thing as the Table 17 i.e. no significant difference between the naviga- 
tion systems when we a\'erage o\ er times of about 20 - 30 minutes. 

In Appendix F there are 6 plots for points 1 - 2, 2 - 1 but not averaged for 
the whole length (just plots for specific times between the points). These plots 
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show, as expected, that the U and V components of the current are different for 
different navigation systems. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



From this thesis we can conclude that: 

• Because of the variability in the LORAN (DISPLAY) results and the need 
for a reasonable measure of statistical strength, it is not recommended that 
LORAN (DISPLAY) be used for giving averaged results of under 20 min- 
utes. Indeed, it is preferable to use LORAN (TD) rather than LORAN 
(DISPLAY) whenever possible. Furthemore, LORAN (DISPLAY) results 
not only show a much greater variability than the LORAN (TD) results, but 
also at a 95% confidence interval consistently fail the statistical hypothesis 
that they come from the same population. 

• MINI RANGER results as expected are very consistent. It is felt that good 
ADCP current estimates could be derived using MINI RANGER observa- 
tions for periods as short as 3 minutes. 

• When averaging for 20 - 30 minutes all the navigation systems give the same 
results. 

• Without averaging we recommended that the MINI RANGER be used 
whenexer possible. If this system is not available, LORAN (TD) should be 
the next choice. 

• The theoritical investigation revealed that when time axeraging oxer 20 
minutes or longer, the choice of the filter to be used for the LORAN system 
xvas not critical. The cruise results appear to support this contention. 

• Bottom Tracking results are systematically biased, but hax^e the same 
standard dexiation as the MINI RANGER results. 

In looking to the future xx’ork on this subject it is recommended that: 

• Further xvork is required to understand and resolxe the Bottom Tracking 
bias problem. 

• If this cruise is done again xx'e suggest checking the LORAN filters by mak- 
ing one long run, dix'iding it into 30 minute segments, and using a different 
filter on each segment. 



46 



APPENDIX A. STATISTICAL ANALYSIS FOR THE LC 408 DATA 



The provided numbers are in nanoseconds . 

( 1 ) TIME COAS TA N T 0 SECONDS (NO FIL TERING ) 



STATION 


N 


MEAN 


STANDARD 


MINIMUM 


MAXIMUM 


STD ERROR 








DEVIATION 


VALUE 


VALUE 


OF MEAN 


WHISKY 


1440 


16312005. 09 


46. 85663572 


16311760. 00 


16312180. 00 


1. 23478077 


X-RAY 


1440 


27511618. 59 


22. 27564796 


27511540. 00 


27511720. 00 


0. 58701487 


YANKEE 


1440 


42735925. 76 


24. 94991412 


42735810. 00 


42736030. 00 


0. 65748797 




(2 ) TIME CONSTANT 5 SECONDS (GOOD RESPONSE ) 




STATION 


N 


MEAN 


STANDARD 


MINIMUM 


MAXIMUM 


STD ERROR 








DEVIATION 


VALUE 


VALUE 


OF MEAN 


WHISKY 


1440 


16312033. 94 


51. 75565604 


16311860. 00 


16312180. 00 


1. 36388129 


X-RAY 


1440 


27511616. 71 


19. 82896815 


27511560. 00 


27511680. 00 


0.52253919 


YANKEE 


1440 


42735928. 43 


23. 39278554 


42735860. 00 


42736020. 00 


0. 61645403 




(3 j TIME CONSTANT 10 SECONDS (SLOWER RESPONSE) 


STATION 


N 


MEAN 


STANDARD 


MINIMUM 


MAXIMUM 


STD ERROR 








DEVIATION 


VALUE 


VALUE 


OF MEAN 


WHISKY 


1440 


16311981. 04 


49. 40674807 


16311790. 00 


16312110. 00 


1. 30198213 


X-RAY 


1440 


27511616. 56 


13. 76653231 


27511580. 00 


27511680. 00 


0. 36277998 


YANKEE 


1440 


42735918. 06 


18. 02041879 


42735850. 00 


42735990. 00 


0. 47487973 




f4) TIME CONSTANT 20 


SECONDS 


(MONITOR 


APPLICATIONS) 












STATION 


N 


MEAN 


STANDARD 


MINIMUM 


MAXIMUM 


STD ERROR 








DEVIATION 


VALUE 


VALUE 


OF MEAN 


WHISKY 


1417 


16311990. 54 


46. 87615596 


16311800. 00 


16312130. 00 


1. 24528015 


X-RAY 


1417 


27511618. 80 


12. 99066592 


27511570. 00 


27511680. 00 


0. 34510122 


YANKEE 


1417 


42735915. 31 


16. 74596778 


42735850. 00 


42735980. 00 


0. 44486202 
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APPENDIX B. TABLES FOR STANDARD DEVIATIONS AND 
COVARIANCES OF VELOCITY COMPONENTS FOR LORAN C (LC 408) 



Table 18. TIME INTERVAL BETWEEN POSITIONS 3 MINUTES 



FILTER 


Gy 


<^v 


<^y,u 


0 sec 


1 1.6 cm; sec 


6.8 cm/sec 


-61.184 cmV sec^ 


5 sec 


10.8 cm. sec 


6.1 cm/ sec 


-51.047 cm^l scc^ 


10 sec 


8.2 cm sec 


4.3 cm/scc 


-27.447 cm^l sec- 


20 sec 


7.7 cm /sec 


4.1 cm/ sec 


-24.029 cm-/ sec^ 



Table 19. TIME INTERVAL BETWEEN POSITIONS 6 MINUTES 



FILTER 


Gy 






0 sec 


5.8 cm sec 


3.4 cm /sec 


-15.296 cm-j sec- 


5 sec 


5.4 cm 'sec 


3.0 cm/sec 


-12.761 cm-j sec- 


10 sec 


4.1 cm sec 


2.1 cm /sec 


-6.861 cm^ j sec^ 
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20 sec 


3.8 cm/sec 


2.0 cm/sec 


-6.007 cm^l sec^ 


Table 20. TIME INTERVAL BETWEEN POSITIONS 12 MINUTES 


FILTER 


Cv 


<^u 


^V, V 


0 sec 


2.9 cm, sec 


1.7 cm/sec 


— 3.824c/? 27 sec^ 


5 sec 


2.7 cm sec 


1 .5 cm/sec 


— 3.190c/?F/ sec" 


10 sec 


2.1 cm sec 


1.1 cm/sec 


— 1.715c/? 7V sec" 


20 sec 


1 .9 cm sec 


1 .0 cm, sec 


— 1 ..502cmV sec’ 


Table 21. TIME INTERVAL BETWEEN POSITIONS 20 MINUTES 


FILTER 


Cv 


Ou 




0 sec 


1.7 cm sec 


1 .0 cm /sec 


— 1.3760/77^ sec- 


5 sec 


1.6 cm,' sec 


0.9 cm, sec 


— 1.148c//77 scc‘ 
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10 sec 


1.2 cm/sec 


0.7 cm/sec 


— 0.6 17cm V sec^ 


20 sec 


1.1 cm/sec 


0.6 cm/sec 


— 0.540cmV sec^ 



r 



Table 22. TIME INTERVAL BETWEEN POSITIONS 40 MINUTES 



FILTER 


Oy 




<^y.u 


0 sec 


0.9 cm sec 


0.5 cm/sec 


— 0.344cmV sec^ 


5 sec 


O.S cm sec 


0.4 cm/sec 


— 0.287cmV sec‘ 


10 sec 


0.6 cm sec 


0.3 cm sec 


— 0. 1 5Acnvj sec^ 


20 sec 


0.5 cm sec 


0.3 cm/sec 


— 0.135ca?F/ sec‘ 
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APPENDIX C. PROGRAM DRIVER FORTRAN 



A. GENERAL REMARKS 

This program calculate the hyperbolic coordinates of a ship expressed in 
LORAN time differences from geographic potisions forward computation , and 
the inverse computation. 

Forward computation; 

• Input parameters - ALAT, AEON, L, DEL, ISYS, IB, IVEL. 

• Output parameters - TD, BETA, BEEN, BLEM, lER. 

Inverse computation: 

• Input parameters - ALAT, AEON, L, DEL, TD. ISYS, ISENT, ISW. IB. 
IVEL. 

• Output parameters - ALAT, ALON. BETA, BEEN, BLEM, ISENT. lER, 
IT. 

B. DESCRIPTION OF PARAMETERS: 

ALAT - Latitude array in seconds of arc: 

• ALAT(l) - Latitude of master station. 

• ALAT(2) - Latitude of slave 1 station. 

• ALAT(3) - Latitude of master station. 

• ALAT(4) - Latitude of slave 2 station. 

• ALAT(5) - Latitude of ship position. 

• ALON(l) - Longitude of master station. 

• ALON(2) - Longitude of slave 1 station. 

• ALON(3) - Longitude of master station. 

• ALON(4) - Longitude of slave 2 station. 

• ALON(5) - Longitude of ship position. 

L - Spheroid code array. 

L(l) - Spheroid code array for first rate: 

• L(l) = 1 WGS 1984. 
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L(l) = 2 


Bessel spheroid. 


L(l) = 3 


Clarke 1858 spheroid. 


L(l) = 4 


Clarke 1866 spheroid. 


L(l) = 5 


Clarke 1880 spheroid. 


L(l) = 6 


Everest spheroid. 


L(l) = 7 


Fischer spheroid. 


L(l) = 8 


International spheroid. 


L(l) = 9 


WGS 1972 spheroid. 


L(2) - Spheroid code for second rate: 


L(2) = 1 


WGS 1984. 


L(2) = 2 


Bessel spheroid. 


L(2) = 3 


Clarke 1858 spheroid. 


L(2) = 4 


Clarke 1866 spheroid. 


L(2) = 5 


Clarke 1880 spheroid. 


L(2) = 6 


E\erest spheroid. 


L(2) = 7 


Fischer spheroid. 


L(2) = 8 


International spheroid. 


L(2) = 9 


WGS 1972 spheroid. 


DEL - Coding delay array, in microseconds: 


DEL(l) 


Coding delay for first rate. 


DEL(2) 


Coding delay for second rate. 



TD - Time difference array, in microseconds: 

• TD(1) Time difference for first rate. 

• TD(2) Time difference for second rate. 



I SYS - Loran system array; 

• ISYS(l) Loran system for first rate. 

1SYS(I ) = A: Loran A 
ISYS(l) =c: Loran C 

• IS'i'S(2) Loran svstem for second rate. 

1SYS(2)'=A: Loran A 
lSYS(2)=c; Loran C 
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BETA - Baseline delay array; 

• BETA(l) Baseline delay for first rale. 

• BETA(2) Baseline array for second rate. 

BEEN - Minimum time difference array, in microseconds: 

• BLEN(l) Minimum time difference for first rate. 

• BLEN(2) Minimum time difference for second rate. 

BLEM - Maximum time difference array, in microseconds: 

• BLEM(l) Maximum time difference for first rate. 

• BLEM(2) Maximum time difference for second rate. 

ISENT - Sentinel for DR position: 



1SENT = 0 


DR position not given. 


ISENT =1 


DR position given. 


ISW Operation code: 


IS\V=1 


Forward computation. 


ISW = 2 


Inverse computation. 


IB - Switch for baseline computation; 


1B=0 


Baseline computation off. 


IB = 1 


Baseline computation on. 


lER - Error code: 


1ER = 1 


In\ alid velocity of propagation code. 


1ER = 2 


Invalid spheroid code. 


IER = 3 


Invalid configuration. 


1ER = 4 


Invalid operation code. 


IER = 5 


Invalid readings. 


IER = 6 


Invalid DR position. 


IT - Iteration code. 



IVEL - Velocity of propagation code array: 

IVEL(l) - Velocity of propagation code for first rate: 
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• 0 = 299.6929 kilometers per microsecond 

• 1 = 299.691 1624 kilometers per microsecond 

• 2 = 299.708 kilometers per microsecond 



IVEL(2) - VElocity of propagation code for second rate: 

• 0 = 299.6929 kilometers per microsecond 

• 1 = 299.6911624 kilometers per microsecond 

• 2 = 299.708 kilometers per microsecond 

* PROGRAM DRIVLR FORTRAN SOURCE CODE * 



C 



10 



20 



DESCRIPTION OF PARAMETERS 
DATA DECLARATION 
IMPLICIT REAL*8(A-H,0-Z) 

REAL*8 ALAT( 5 ) , AL0N( 5 ) , DEL( 2 ) ,TD( 2) , BETA( 2 ) , BLEN( 2 ) , BLEM( 2) , LN 
CHARACTER*! ISYS(2) 

INTEGER H,M,L(2),IVEL(2),ISYS(2) 

ALAT( 1)=3600. 0D0*39. ODO + 33. 0D0*60. ODO + 

+ 46. 0D0*60. ODO + 

+ 33. 0D0*60. ODO + 

+ 19. 0D0*60. ODO + 

+ 45. 0D0*60. ODO + 



ALAT(2)=3600. 0D0*38. ODO 
ALAT(3)=3600. 0D0*39. ODO 
ALAT(4)=3600. 0D0*35. ODO 
ALAT(5)=3600. 0D0*36. ODO 
AL0N(1)=3600. 0D0*118. ODO 
ALON(2)=3600. 0D0*122. ODO 
AL0N(3)=3600. 0D0*118. ODO 
AL0N(4)=3600. 0D0*114. ODO 
ALON(5)=3600. 0D0*121. ODO 
L(l) = 1 
L(2) = 1 

DEL(l) = 27000. ODO 
DEL(2) = 40000. ODO 
ISYS( 1)=1 
ISYS(2)=1 
ISENT=1 



+ 49. 0D0*60. ODO + 
+ 29. 0D0*60. ODO + 
+ 49. 0D0*60. ODO + 
+ 48. 0D0*60. ODO + 
+ 55. 0D0*60. ODO + 



6. 740D0 
57. IIODO 
6. 740D0 
18. 3050D0 
0. ODO 
55. 8160D0 
43. 9750D0 
55. 8160D0 
16. 8810D0 
0. ODO 



ISW=2 

IVEL(1)=1 

IVEL(2)=1 

IB=1 



CALL LORAN ( ALAT , ALON , L , DEL , TD , I SYS , BETA , BLEN , BLEM , I SENT , I SW , I B , 
*IER,IT,IVEL) 

DO 200 1=1,914 

READ(7,10) H,M,S,TD(1),TD(2) 

F0RMAT(I2,1X,I2,1X,F5. 2,21X,F8. 2,1X,F8. 2) 

IER=0 

CALL LORAN ( ALAT , ALON , L , DE L , TD , I S YS , BETA , BLEN , BLEM , I SENT , I SW , I B , 
*IER,IT,IVEL) 

WRITE(9,20) I,H,M,S,ALAT(5),ALON(5),IT,IER 
F0RMAT(I3,1X,I2,1X,I2,1X,F5. 2,2X,F10. 2,2X,F10. 2 , 2X, 12 , 2X, I 1) 
F=ALAT(5)/3600. ODO 
LN=-ALON(5)/3600. ODO 
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WRITE(10,30) I,H,M,S,F,LN 

30 F0RMAT(I3,1X,I2,1X,I2,1X,F5. 2,2X,F10. 7,2X,F12. 7) 

READ(8,40) F1,LNI 
40 FORMAT(9X,F10. 7,1X,F12. 7) 

DF=F-F1 

DL=LN-LN1 

WRITE(6,*) DF,DL,LN1 
200 CONTINUE 
STOP 
END 

SUBROUTINE LORAN ( ALAT,AL0N, L, DEL, TD, ISYS , BETA, BEEN, BLEM, 
*ISENT,ISW,IB,IER,IT,IVEL) 

C 

C DATA DECLARATION 

C 

IMPLICIT REAL*8(A-H,0-Z) 

INTEGER SOUTH, EAST, WEST 
C 

C STORAGE ALLOCATION 

C 

DIMENSION DEL(2) ,TD(2) ,BETA(2) ,BLEN(2) ,BLEM(2) ,ISYS(2) ,H(2) ,R(2) , 
*IVEL(2),VELL(2),L(2) 

DIMENSION VEL(3) 

DIMENSION TEMP1(4),TEMP2(4) 

DIMENSION ALAT(5) ,ALON(5) 

DIMENSION A(9),B(9) 

DATA A/ 6378137. ODO , 6377397. 1550D0 , 6378293. 6450D0 , 

* 6378206. 4000D0 , 6378249. 1450D0, 6377276. 3450D0, 6378166. OOOODO, 
*6378388. OOOODO , 6378135. OOOODO/ 

DATA B/6356752. 31420D0 ,6356078. 9628D0, 6356617. 9376D0, 

*6356583. 8000D0, 63565 14. 8695D0 , 6356075. 4131D0 , 6356784. 2836D0, 
*6356911. 946100,6356750. 5200D0/ 

DATA Al/24. 0305D0/ ,A2/-0. 40758D0/ , A3/0. 00346776D0/ 

DATA Bl/0. 510483D0/,B2/-0. 011402D0/ ,B3/0. 001760D0/ 

DATA ARCl/0. 484813681110D-5/ 

DATA LORANC/1/ 

DATA FAZ/0. ODO/ , IK/0/ ,TEMPl/4*0. ODO/ ,TEMP2/4*0. ODO/ 

DATA ZERO/0. ODO/ 

DATA NORTH , SOUTH , EAST, WEST/ IHN , IHS , IHE , IHW/ 

DATA MINUS/ IH-/ 

DATA VEL/299. 692900,299. 691162400,299. 708D0/ 

DATA BAZ/0. 0D0/,JER/0/ 

DATA ALATS/0. ODO/ ,ALONS/0. 0D0/,N/0/ 

C 

C INITIALIZE ERROR CODE 

C 

IER=0 

C 

C TEST SWITCH FOR BASELINE COMPUTATION 

C 

IF (IB) 5,150,5 
C 

C BASELINE COMPUTATION 

C 

C 
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TEST FOR VALID VFLOCITY OF PROPAGATION CODE FOR BOTH CHAINS 

5 DO 10 1=1,2 

IF (IVEL(I).GT. 2) GO TO 15 
10 CONTINUE 
GO TO 20 
15 IER=1 

GO TO 320 

DETERMINE VELOCITY OF PROPAGATION FOR BOTH CHAINS 

20 DO 25 1=1,2 

IVEL(I)=IVEL(I)+1 

KK=IVEL(I) 

25 VELL(I)=VEL(KK) 

A4A=VELL(1) 

A4B=VELL(2) 

I, _ 

TEST FOR CHANGE IN SPHEROID CODE, FOR VALID SPHEROID' (JODE 
AND DETERMINE APPROPRIATE SPHEROID COSTANTS 

DO 30 J=l,2 

IF (L(J).GE. l.AND. L(J).LE. 9) GO TO 30 
GO TO 55 
30 CONTINUE 

IF (L(l).EQ. L(2)) GO TO 35 
GO TO 55 
35 K=L(1) 

IF (IK-K) 40,60,40 
40 DO 50 1=1,9 

IF (I-K) 50,45,50 
45 A5=A(I) 

A6=B(I) 

A55=A5*A5 

A10=(A55-A6*A6)/A55 

IK=K 

GO TO 60 
50 CONTINUE 
55 IER=2 

GO TO 320 

TEST FOR CHANGE IN FIXED STATIONS 
60 DO 80 1=1,4 

IF (TEMPl(I)-ALAT(I)) 70,65,70 
65 IF (TEMP2(I)-ALON(I)) 70,80,70 
70 DO 75 J=l,4 

TEMP1(J)=ALAT(J) 

75 TEMP2(J)=ALON(J) 

GO TO 85 
80 CONTINUE 
GO TO 145 

TEST FOR TRIAD OR TETRAD CONFIGURATION AND ACTIVATE 
APPROPRIATE SENTINEL 
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85 A12=ALAT( 1) -ALAT( 3)+AL0N( 1) -ALON( 3) 

A12=DABS(A12) 

A13=ALAT(2)-ALAT(4)+AL0N(2)-AL0N(4) 

A13=DABS(A13) 

IF (A12-0. OOIDO) 90,90,105 
90 IF (A13-0. OOIDO) 95,95,100 
95 IER=3 

GO TO 320 
100 All=-1. ODO 
GO TO 110 
105 All=l. ODO 

TEST FOR VALID OPERATION CODE 

110 IF (ISW-1) 120,125,115 
115 IF (ISW-2) 125,125,120 
120 IER=4 

GO TO 320 

COMPUTE GEODETIC DISTANCE AND SALT-WATER RETARDATION FROM 
MASTER 1 STATION TO SLAVE 1 STATION 

125 CALL SODIN (ALAT( 1) ,ALON( 1) ,ALAT(2) ,AL0N(2) ,L,A35 ,FAZ,BAZ, JER) 
AA10=A35/A4A 
M=1 
1=1 

GO TO 295 
130 AA11=A47 

COMPUTE GEODETIC DISTANCE AND SALT-WATER RETARDATION FROM 
MASTER 2 STATION TO SLAVE 2 STATION 

CALL SODIN (ALAT(3),AL0N(3),ALAT(4),AL0N(4),L,A35,FAZ,BAZ,JER) 
CC10=A35/A4B 
M=2 
1=2 

GO TO 295 
135 CC11=A47 

COMPUTATION OF BASELINE DATA 
BETA(1)=AA10+AA11 
BETA(2)=CC10+CC11 
DO 140 1=1,2 
BLEN(I)=DEL(I) 

140 BLEM(I)=2. 0D0*BETA(I)+DEL(I) 

TURN SWITCH FOR BASELINE COMPUTATION OFF 

145 IB=0 
RETURN 

TEST SWITCH FOR FORWARD OR INVERSE COMPUTATION 
150 GO TO (155,190),ISW 

FORWARD COMPUTATION 
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COMPUTE GEODETIC DISTANCE AND SALT-WATER RETARDATION FROM 
SHIP TO MASTER 1 STATION 

155 CALL SODIN ( ALAT(5 ) ,AL0N(5 ) , ALAT( 1) ,ALON( 1) ,L,A35 ,FAZ ,BAZ , JER) 
A76=A35/A4A 
A72=A76 
M=1 
1=7 

GO TO 295 
160 A77=A47 
A73=A77 

COMPUTE GEODETIC DISTANCE AND SALT-WATER RETARDATION FROM 
SHIP TO SLAVE 1 STATION 

CALL SODIN (ALAT(5),AL0N(5),ALAT(2),AL0N(2),L,A35,FAZ,BAZ,JER) 
A74=A35/A4A 
M=1 
1=8 

GO TO 295 
165 A75=A47 

COMPUTE GEODETIC DISTANCE AND SALT-WATER RETARDATION FROM 
SHIP TO SLAVE 2 STATION 

CALL SODIN (ALAT(5),AL0N(5),ALAT(4),AL0N(4),L,A35,FAZ,BAZ,JER) 

A70=A35/A4B 

M=2 

1=9 

GO TO 295 
170 A71=A47 

COMPUTE GEODETIC AND SALT-WATER RETARDATION FROM 
SHIP TO MASTER 2 STATION 

IF (All) 185,175,175 

175 CALL SODIN ( ALAT(5 ) ,ALON(5 ) , ALAT( 3) , AL0N( 3) ,L, A35 ,FAZ , BAZ , JER) 
A72=A35/A4B 
M=2 
1=10 

GO TO 295 
180 A73=A47 

COMPUTE LORAN TIME DIFFERENCES OF FIX 

185 TD(1)=AA10+AA11+A74+A75-A76-A77+DEL(1) 
TD(2)=CC10+CC11+A70+A71-A72-A73+DEL(2) 

RETURN 

INVERSE COMPUTATION 



STORE DR POSITION OF SHIP 
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190 ALATDR=ALAT( 5 ) 

AL0NDR=AL0N(5) 

TEST FOR VALID MICROSECOND READINGS 
DO 195 1=1,2 

IF (TD(I). LE. BLEM(I). AND. TD(I). GE. BLEN(I)) GO TO 195 
IER=5 
GO TO 320 
195 CONTINUE 

CONVERT SHIPS POSITION TO RADIANS 

IF (ISENT) 205,205,200 
200 ALATR=ALAT(5)*ARC1 
AL0NR=AL0N(5)*ARC1 

INITIALIZE lETRATION COUNTER 

205 IT=0 

GO TO 215 

TEST ITERATION COUNTER 

210 IF (IT. LT. 30) GO TO 215 
IER=6 
GO TO 320 

COMPUTE SINES AND COSINES OF SHIPS POSITION 

215 SINS=DSIN(ALATR) 

COSS=DCOS(ALATR) 

COMPUTE GEODETIC DISTANCE, SIN AND COSINE OF AZIMUTH, AND 
SALT-WATER RETARDATION FROM SHIP TO SLAVE 2 STATION 

CALL SODIN (ALAT(5),ALON(5),ALAT(4),ALON(4),L,A35,FAZ,BAZ,JER) 

M=2 

1=3 

GO TO 280 
220 C1=A35 
D1=C1/A4B 
C2=DSIN(FAZ) 

C3=DCOS(FAZ) 

C101=A47 

COMPUTE GEODETIC DISTANCE, SIN AND COSINE OF AZIMUTH, AND 
SALT-WATER RETARDATION FROM SHIP TO MASTER 2 STATION 

CALL SODIN (ALAT(5),ALON(5),ALAT(3),ALON(3),L,A35,FAZ,BAZ,JER) 

M=2 

1=4 

GO TO 280 
225 C4=A35 
D4=C4/A4B 
C5=DSIN(FAZ) 
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C6=DC0S(FAZ) 

C104=A47 

COMPUTE GEODETIC DISTANCE, SIN AND COSINE OF AZIMUTH, AND 
SALT-WATER RETARDATION FROM SHIP TO SLAVE 1 STATION 

CALL SODIN (ALAT(5),AL0N(5),ALAT(2),AL0N(2),L,A35,FAZ,BAZ,JER) 

M=1 

1=5 

GO TO 280 
230 C7=A35 
E7=C7/A4A 
C8=DSIN(FAZ) 

C9=DC0S(FAZ) 

C107=A47 

IF (All) 245,235,235 

COMPUTE GEODETIC DISTANCE, SIN AND COSINE OF AZIMUTH, AND 
SALT-WATER RETARDATION FROM SHIP TO MASTER 1 STATION 
(FOR TETRAD CONFIGURATIONS ONLY) 

235 C10=C7 
C11=C8 
C12=C9 

CALL SODIN (ALAT(5),ALON(5),ALAT(l),ALON(l),L,A35,FAZ,BAZ,JER) 
M=1 
1=6 

GO TO 280 
240 C7=A35 
D7=C7/A4A 
C8=DSIN(FAZ) 

C9=DCOS(FAZ) 

C108=A47 

DETERMINATION OF TIME DIFFERENCE COSTANTS AND DIFFERENTIAL 

CORRECTIONS IN Y AND X 

(FOR TETRAD CINFIGURATIONS ONLY) 

C13=TD(2)-CC10-CC11-C101+C104-DEL(2) 

C17=C13*A4B 

C18=TD(1)-AA10-AA11-C107+C108-DEL(1) 

C22=C18*A4A 

C23=C1-C17 

C24=C4 

C25=C7+C22 

C26=C10 

C27=((C2-C5)*(C25-C26)*(C23-C24)*(C11-C8)) 

C29=((C2-C5)*(C12-C9)*(C3-C6)*(C8-C11)) 

C30=C27/C29 

C28=(C23-C24+C30*(C3-C6))/(C5-C2) 

GO TO 250 

DETERMINATION OF TIME DIFFERENCE COSTANTS AND DIFFERENTIAL 

CORRECTIONS IN Y AND X 

(FOR TRIAD CINFIGURATIONS ONLY) 
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245 C13=TD(2)-CC10-CC11-C101+C104-DEL(2) 

C17=C13*A4B 

C18=TD( 1)-AA10-AA11-C107+C104-DEL( 1) 

C22=C18*A4A 

C23=C1-C17 

C24=C4 

C25=C7-C22 

C27=C2*(C25-C24)+C23*(C5-C8)+C8*C24-C5*C25 

C29=C2*(C6-C9)+C3*(C8-C5)+C5*C9-C8*C6 

C30=C27/C29 

C28=(C23-C24+C30*(C3-C6))/(C5-C2) 

DETRMINATION OF DIFFERENTIAL CORRECTIONS IN LATITUDE AND 
LONGITUDE AND PRELIMINARY GEOGRAPHIC POSITION OF SHIP 

250 C100=l. 0D0-A10*SINS*SINS 
C102=DSQRT(C100) 

C32=A5/C102 

C31=C32*(( 1. 0D0-A10)/C100) 

C33=C30/C31 

C34=(-C28/(C32*C0SS)) 

ALATR=ALATR+C33 

AL0NR=AL0NR+C34 

ALAT(5)=ALATR/ARC1 

AL0N(5)=AL0NR/ARC1 

WRITE(13,*) ALAT(5),AL0N(5),'ALAT ALON FM LORAN' 

TEST OF CONVERGENCE OF SOLUTION 
IF (All) 260,255,255 

255 H(1)=AA10+AA11+E7+C107-D7-C108+DEL(1) 

GO TO 265 

260 H( 1)=AA10+AA11+E7+C107-D4-C104+DEL( 1) 

265 H(2)=CC10+CC11+D1+C101-D4-C104+DEL(2) 

DO 270 1=1,2 
270 R(I)=TD(I)-H(I) 

ANUM=DABS(R( 1) )+DABS(R(2) ) 

IF (ANUM. LT. 0. OOOIDO) GO TO 275 
IT=IT+1 
GO TO 210 
275 RETURN 

CONVERSION OF AZIMUTHS FROM NORTH TO SOUTH ORIENTATION 

280 IF (FAZ-648000. ODO) 285,285,290 
285 FAZ=(FAZ+648000. 0D0)*ARC1 
GO TO 295 

290 FAZ=( FAZ-648000. 0D0)*ARC1 



FORWARD AND INVERSE COMPUTATION 
OF SALT-WATER RETARDATION (FOR LORAN C ONLY) 

295 IF (ISYS(M). EQ. LORANC) GO TO 300 
A47=0. ODO 
GO TO 315 
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300 A46=A35/1609. 349D0 

IF (A46-100. ODO) 305,310,310 
305 A47=(B1/A46)+B2+B3*A46 
GO TO 315 

310 A47=(A1/A46)+A2+A3*A46 

BRANCHING FOR QUANTITIES TO BE USED IN SOLUTION 

315 GO TO (130,135,220,225,230,240,160,165,170,180),! 

PRINT ROUTINE FOR INVALID DATA 

320 IF (lER. LE.4) GO TO 340 

CALL ANGLE (N,ALONDR,ISGN,LOND,LONM,ALONS) 

CALL ANGLE (N,ALATDR,ISGM,LATD,LATM,ALATS) 

IF (ISGM.EQ. MINUS) GO TO 325 
KLA=N0RTH 
GO TO 330 
325 KLA=S0UTH 

330 IF (ISGN.EQ. MINUS) GO TO 335 
KL0=WEST 
GO TO 340 
335 KLO^EAST 

340 GO TO (345, 350, 355, 360, 365, 370), lER 
345 PRINT 375 
RETURN 

350 PRINT 380 
RETURN 

355 PRINT 385 
RETURN 

360 PRINT 390 
RETURN 

365 PRINT 395,LATD,LATM,KLA,LOND,LONM,KLO,TD(l) ,TD(2) 

RETURN 

370 PRINT 400,LATD,LATM,KLA,L0ND,LONM,KL0,TD(l),TD(2) 

RETURN 

FORMAT STATEMENTS 

375 FORMATdHl, IX, 'INVALID VELOCITY OF PROPAGATION CODE, DATA NOT PROC 
*ESSED' ) 

380 FORMATdHl, IX, 'INVALID SPHEROID CODE, DATA NOT PROCESSED') 

385 FORMAT( IHl, IX, ' INVALID CONFIGURATION, DATA NOT PROCESSED') 

390 FORMATdHl, IX, 'INVALID OPERATION CODE, DATA NOT PROCESSED') 

395 F0RMAT(18X,2I3,A1,6X,I4,I3,A1,16X,F9. 2,13X,F9. 2,3X, ' INVALID READIN 
*GS ' ) 

400 F0RMAT(18X,2I3,A1,6X,I4,I3,A1,16X,F9. 2,13X,F9, 2,3X, 'INVALID DR POS 
*ITI0N' ) 

END 



SUBROUTINE SODIN ( ALATl ,AL0N1 ,ALAT2 , ALON2 ,L,S ,FAZ ,BAZ , lER) 
C 

IMPLICIT REAL*8(A-H,0-Z) 

C 

DIMENSION X(2),E(2),F(2),L(2) 
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DATA PI/3. 14159265359D0/,ARCl/0. 484813681110D-5/ 
DATA A,B/6378137. ODO, 6356752. 3142D0/ 

INITIALIZE ERROR CODE AND DETERMINE COSTANTS 

A14=l. 0D0-(B/A) 

A56=A14*A14 
A50=l. 0D0+A14+A56 
A51=A50-1. ODO 
A52=A56/2. ODO 
A53=-A51/2. ODO 
A54=A56/16. ODO 
A55=A56/8. ODO 
A57=A56*l. 250D0 
A58=A56/4. ODO 

TEST FOR IDENTICAL STATIONS 

30 IF(ALAT1-ALAT2) 45,35,45 
35 IF(ALONl-ALON2) 45,40,45 
40 IER=2 

GO TO 255 

COMPUTATION OF PARAMETRIC LATITUDE OF STATIONS 

45 X1=ALAT1*ARC1 
Y1=AL0N1*ARC1 
X2=ALAT2*ARC1 
Y2=AL0N2*ARC1 
X(1)=X1 
X(2)=X2 
DO 50 1=1,2 

BETA=DATAN( ( 1 . ODO -A 14 ) *DS IN( X( I ) ) /DCOS ( X( I ) ) ) 
E(I)=DSIN(BETA) 

50 F(I)=DCOS(BETA) 

COMPUTATION OF SPHERICAL DIFFERENCE OF LONGITUDES 

A59=-Y2 

A60=-Y1 

C35=A59-A60 

C36=DABS(C35) 

IF(C36-PI) 60,55,55 
55 A16=2. 0D0-'-PI-C36 

GO TO 65 
60 A16=C36 

65 IF(A16) 75,70,75 

70 A16=0.5D-7 

COMPUTATION OF GEODETIC DISTANCE 

75 A17=DSIN(A16) 

A18=DCOS(A16) 

A19=E(1)*E(2) 

A20=F(1)*F(2) 

A21=A19+A20*A18 
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A80=A17*F(2) 

A81=A80*A80 

A82=E(2)*F(1)-E(1)*F(2)*A18 

A83=A82*A82 

A84=A81+A83 

A22=DSQRT(A84) 

A23=(A20*A17)/A22 
A24=l. 0D0-A23*A23 
A90=A22/A21 
A25=DATAN(A90) 

IF(A25) 80,85,85 
'■ 80 A25=(PI/2. 0D0)+A25 

85 A26=A25*A25 

A27=l. 0D0/A22 
A91=A21/A22 
A29=A24*A24 

A30=( A50*A25 )+A 19*( A5 1*A22-A52*A26*A27 ) 
A31=A24*(A53*A25+A53*A22*A21+A52*A26*A91) 

A3 2=A 19*A 1 9*( - A5 2*A2 1*A2 2 ) 

A33=A29*(A54*A25+A54*A22*A21-A52*A26*A91-A55*A22*A21**3) 
A34=A19*A24*(A52*A26*A27+A52*A22*A21*A21) 

S=( A30+A3 1+A32+A33+A34 )*B 

COMPUTATION OF DIFFERENCE OF LONGITUDE ON THE REDUCED SPHERE 

A3 6=( A5 1*A25 ) +A 1 9 * ( - A5 2*A2 2 --A5 6*A 2 6*A2 7 ) 
A37=A24-'-(-A57*A25+A58*A22*A21+A56*A26*A91) 
A38=(A36+A37)*A23+A16 

COMPUTATION OF GEODETIC FORWARD AZIMUTH 

A39=DSIN(A38) 

A40=DC0S(A38) 

A4 1=( E( 2 )*F( 1 ) -A40*E( 1 )*F( 2 ) ) / ( A39*F( 2 ) ) 

IF(A41) 95,90,95 
90 A41=0.5D0-7 

95 A42=l. 0D0/A41 

A43=DATAN(A42) 

IF(C35) 120,100,100 
100 IF(C35-PI) 105,115,115 
105 IF(A41) 110,140,140 
110 A43=PI+A43 
GO TO 140 

115 IF(A41) 130,135,135 
120 IF(C35+PI) 105,105,125 
125 IF(A41) 130,135,135 
130 A43=PI-A43 
GO TO 140 

135 A43=2. 0D0*PI-A43 
140 A43=A43+PI 

A43=A43-2. 0D0*PI 
IF(A43) 145,150,150 
145 A43=A43+2. 0D0*PI 
150 FAZ=A43/ARC1 

FAZ=FAZ+648000. ODO 
IF(FAZ-1296000. ODO) 160,155,155 
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155 FAZ=FAZ-1296000. ODO 

COMPUTATION OF BACK AZIMUTH NOT INCLUDED IN THIS VERSION 

160 BAZ=0.0D0 
RETURN 

PRINT ROUTINE FOR INVALID DATA 
255 PRINT 270 

CALL ANGLE ( 3 , ALAT 1 , I SGN 1 , LAT ID , LAT IM , ALAT 1 S ) 

CALL ANGLE ( 3 , ALON 1 , I SGN2 , LON ID , LON IM , ALON 1 S ) 

CALL ANGLE ( 3 , ALAT2 , ISGN3 , LAT2D , LAT2M , ALAT2S ) 

CALL ANGLE ( 3 , AL0N2 , I SGN4 , L0N2D , L0N2M , AL0N2S ) 

265 PRINT 280,ISGN1,LAT1D,LAT1M,ALAT1S,ISGN2,LON1D,LON1M,ALON1S,ISGN3, 
*LAT2D , LAT2M , ALAT2S , I SGN4 , L0N2D , L0N2M , AL0N2S 
RETURN 

FORMAT STATEMENTS 

270 FORMATdHl, 5 IX, 'SUBROUTINE SODIN INVALID DATA'/) 

280 FORMATC IX, 4(2X,A1, 213, F7. 3), 2X, 12, 2X, 'IDENTICAL STATIONS') 

END 



SUBROUTINE ANGLE (N,ARGS , ISIGN, IDEG,MIN,SEC) 



IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION ASEC(6) 

DATA ASEC/59. 45D0,59. 945D0,59. 9945D0,59. 99945D0,59. 999945D0, 
*59. 9999945D0/ 

DATA MINUS/ 1H-/,IBLANK/1H / 

ANGLE CONVERSION 

IF (ARGS) 10,5,5 
5 ISIGN=IBLANK 
BRGS=ARGS 
GO TO 15 
10 ISIGN=MINUS 
BRGS=-ARGS 

15 IDEG=BRGS/3600. ODO 
ARGT=IDEG*3600 
AMIN=BRGS-ARGT 
MIN=AMIN/60. ODO 
AMINT=MIN*60 
SEC=AMIN-AMINT 
M=N+1 

IF (SEC-ASEC(M)) 30,20,20 
20 SEC=0. ODO 
MIN=MIN+1 

IF (MIN-60) 30,25,25 
25 MIN=0 

IDEG=IDEG+1 
30 RETURN 
END 
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APPENDIX D. TABLE FOR STANDARD DEVIATIONS AND 
COVARIANCES OF VELOCITY COMPONENTS FOR MINI RANGER 



TIME IN- 
TERVAL 


Gy 




(^v.u 


3 MIN 


1.3 cm/sec 


1.3 cm/sec 


0.1234 cnv/ sec^ 


6 MIN 


0.64 cm sec 


0.64 cm/sec 


0.0386 cm^/ sec^ 


12 MIN 


0.32 cm, sec 


0.32 cm,/ sec 


0.01 16 c//?V sec^ 


20 MIN 


0.19 cm sec 


0.19 cm/sec 


10“* X 0.0486c//? V sec^ 


40 MIN 


10"^ X 0 . 98 c / 7 z / sec 


10”^ X 0.97c/??/ sec 


10"^ X 0.l7c//?V sec‘ 


60 MIN 


10“^ X 0.67cm/ sec 


10“^ X 0.66cm/ sec 


10"^ X 0.1080c/??V sec^ 
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APPENDIX E. AVERAGE CURRENT PROFILES 




Figure 17. AVERAGE CURRENT PROFILE FROM POINT 1 TO POINT 3 
USING NAVIGATION DATA FROM MINI RANGER 
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UELOCITY <CM/s) 
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Figure 18. AVERAGE CURRENT PROFILE FROM POINT 3 TO POINT 2 
USING NAVIGATION DATA FROM MINI RANGER 
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OERTH < at> > 



UELOCITY <cM/s) 
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Figure 19. 



24 -18 -12 -6 0 6 12 18 24 30 

-I — ■ — I — ■ — I — I — I 1 I I I 1 I I I I I ■ I 

1 

\ 



.. 

X* 




AVERAGE CURRENT PROFILE FROM POINT 2 TO POINT 3 
USING NAVIGATION DATA FROM MINI RANGER 
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UELOCITY (cm/s) 




Figure 20. AVERAGE CURRENT PROFILE FROM POINT 3 TO POINT 4 
USING NAVIGATION DATA FROM MINI RANGER 
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UELOCITY <cri/s) 
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Figure 21. AVERAGE CURRENT PROFILE FROM POINT 1 TO POINT 3 
USING NAVIGATION DATA FROM LORAN (TD) 
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UELOCITY <cn/s) 



-30 -24 -18 -12 -6 0 6 12 18 24 30 




Figure 22. AVERAGE CURRENT PROFILE FROM POINT 3 TO POINT 2 
USING NAVIGATION DATA FROM LORAN (TD) 
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UELOCITY <cn/s) 



( 
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Figure 23. AVERAGE CURRENT PROFILE FROM POINT 2 TO POINT 3 
USING NAVIGATION DATA FROM LORAN (TD) 
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UELOCITY (cM/s) 
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Figure 24. AVERAGE CURRENT PROFILE FROM POINT 3 TO POINT 4 
USING NAVIGATION DATA FROM LORAN (TD) 
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UELOCITY <c«/s) 



-30 -24 -18 -12 -6 0 6 .12 18 24 30 




Figure 25. AVERAGE CURRENT PROFILE FROM POINT 1 TO POINT 3 
USING NAVIGATION DATA FROM LORAN (DISPLAY^) 
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UELOCITV <c«/s) 
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Figure 26. AVERAGE CURRENT PROFILE FROM POINT 3 TO POINT 2 
USING NAVIGATION DATA FROM LORAN (DISPLAY) 
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UELOCITY <cti/s) 




Figure 27. AVERAGE CURRENT PROFILE FROM POINT 2 TO POINT 3 
USING NAVIGATION DATA FROM LORAN (DISPLAVO 
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UELOCITY <cfi/s) 




Figure 28. AVERAGE CURRENT PROFILE FROM POINT 2 TO POINT 4 
USING NAVIGATION DATA FROM LORAN (DISPLAY) 
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APPENDIX F. CURRENT PROFILES FOR SPECIFIC TIMES 



UELOCITY <ctt/s) 




Figure 29. CURRENT PROFILE FROM POINT 1 TO POINT 2 USING 
NAVIGATION DATA FROM MINI RANGER FOR A SPECIFIC 
TIME (14 33 11) 
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UELOCITY <Ctt/5) 
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Figure 30. CURRENT PROFILE FROM POINT 2 TO POINT 1 USING 
NAVIGATION DATA FROM MINI RANGER FOR A SPECIFIC 
TIME (15 06 10) 
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UELOCITY <cfi/s) 
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Figure 31. CURRENT PROFILE FROM POINT 1 TO POINT 2 USING 
NAVIGATION DATA FROM LORAN (TD) FOR A SPECIFIC 
TIME (14 33 11) 
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UELOCITY <CM/s) 



-30 -24 -18 -12 -6 0 6 12 18 24 30 




Figure 32. CURRENT PROFILE FROM POINT 2 TO POINT 1 USING 
NAVIGATION DATA FROM LORAN (TD) FOR A SPECIFIC 
TIME (15 06 10) 
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UELOCITY <cii/s) 




Figure 33. CURRENT PROFILE FROM POINT 1 TO POINT 2 USING 
NAVIGATION DATA FROM LORAN (DISPLAVO FOR A SPE- 
CIFIC TIME (14 33 II) 
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UELOCITY (cii/s) 



-30 -24 -18 -12 -6 0 6 12 18 24 30 




Figure 34. CURRENT PROFILE FROM POINT 2 TO POINT 1 USING 
NAVIGATION DATA FROM LORAN (DISPLAY) FOR A SPE- 
CIFIC TIME (15 06 10) 
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